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ABSTRACT 


The  transfer  matrix  method  employed  for  vibration  analysis 
of  periodic  structure  in  the  present  study  is  somewhat 
similar  in  principle  to  the  classical  transfer  matrix  method 
which  has  long  been  used  for  a  continuous  beam  on  multiple 
supports  (either  spring  or  knife-edge  supports).  However, 
improvements  are  made  on  the  method  by  combining  the 
classical  transfer  matrix  method  with  the  concept  of  the 
propagation  constant  and  the  Finite  Element  formulation  of 
an  impedance  matrix.  The  Finite  Element  formulation  enables 
one  to  do  the  vibration  analysis  no  matter  how  complex  the 
structure  may  be;  while  the  use  of  transfer  matrices  reduces 
the  computational  effort  virtually  to  that  of  a  single 
periodic  unit. 

The  propagation  constant  which  is  obtained  in  the  form 
of  a  logarithmic  eigenvalue  of  the  transfer  matrix  helps  in 
the  understanding  of  the  vibration  by  indicating  which 
flexural  waves  attenuate  and  which  propagate  along  the 
periodic  structure.  Consequently,  vibration  problems  of 
periodic  structures  which  could  be  very  difficult  to  analyze 
by  other  methods  can  be  more  easily  and  efficiently  analyzed 
by  this  improved  transfer  matrix  method. 

The  numerical  examples  taken  in  the  present  study  are 
mostly  for  uniform  Timoshenko  beams  with  different 
intermediate  supports  and  extreme  boundary  conditions.  The 
natural  frequencies  or  the  steady-state  forced  harmonic 
response  of  such  beams  are  obtained  and  the  effect  of  rotary 
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inertia  and  shear  deformation  on  the  flexural  motions  of  a 
multi-span  beam  is  investigated.  The  application  of  the 
present  method  is  by  no  means  restricted  to  the  uniform  beam 
vibration  problems.  However,  the  numerical  examples  and 
their  results  may  sufficiently  demonstrate  the  versatility 
and  superiority  of  the  present  method  over  other  available 
methods  for  analysis  of  vibration  of  periodic  structures. 
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I.  INTRODUCTION  AND  LITERATURE  SURVEY 


A.  Periodic  Structure 

Many  practical  structures  are  constructed  in  the  form 
of  a  "periodic  structure"  mainly  for  appearance  and  lower 
building  costs.  The  interior  of  a  periodic  structure 
consists  of  identical  substructures  connected  in  identical 
ways.  Figure  1.1  shows  two  simplified  examples  of  periodic 
structures.  In  the  first  example  of  a  simple  periodic 
structure,  the  identical  periodic  unit  simply  repeats  over 
and  over  again,  while  the  periodic  unit  in  the  second 
example  is  itself  a  periodic  structure  composed  of  periodic 
subunits.  Being  the  generalization  of  the  simple  periodic 
structure,  the  latter  is  called  a  compound  periodic 
structure.  The  constructions  of  bridges,  ships  and 
airplanes  contain  numerous  examples  of  the  second  type, 
while  railroad  tracks,  shafts  with  bearings  and  discs  can  be 
considered  good  examples  of  the  first  type. 

For  vibration  analysis,  these  beam  and  plate-like 
structures  with  regularly  spaced  supports  or  stiffeners  are 
often  approximated  by  a  uniform  beam  model  resting  on 
equi-spaced  supports,  known  as  a  'periodic  beam'.  These 
supports  provide  translational  and  rotational  restraints  to 
the  flexural  motion  and  thus  may  be  represented  by 
knife-edges  or  by  translational  and  rotational  elastic 
springs . 
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Being  a  mono-coupled  system,  the  uniform  beam  resting 
on  uniformly  spaced  knife-edge  supports  is  one  of  the 
simplest  case  of  the  periodic  structure,  and  hence  it  has 
been  the  favorite  subject  of  earlier  studies  by  Miles[16], 
Ayre  and  Jacobsen  [1]  and  many  others.  Using  a  graphical 
technique  of  building  up  a  nomograph  of  natural  frequencies 
vs  number  of  spans,  it  was  shown  in  [1]  that  the  natural 
frequencies  of  the  multi-span  beam  can  be  obtained  for  two 
types  of  extreme  boundaries,  namely,  simple  and  clamped. 

The  repetitive  nature  of  the  groups  of  natural  frequencies, 
as  many  frequencies  in  each  group  as  there  are  spans,  was 
observed  in  his  study.  A  difference  equation  method  was 
also  shown  to  be  applicable  to  find  the  natural  frequencies 
of  such  beams  [16].  However,  this  uniform  beam  model 
resting  on  knife-edge  supports  can  not  satisfactorily 
approximate  many  complicated  types  of  periodically  stiffened 
structures  commonly  used  in  aeronautical  and  naval 
f  ramewor ks . 

B.  Transfer  Matrix 

When  the  continuous  beam  is  on  uniformly  spaced  elastic 
springs  which  provide  for  constraints  to  the  deflection  and 
rotation  of  the  beam  the  system  becomes  mult i -coupled  and 
requires  matrix  representation  of  coupling  displacements  and 
forces.  For  this  type  of  problem  the  transfer  matrix  which 
relates  such  coupling  displacements  and  forces  in  successive 
substructures  was  efficiently  used  by  Lin[10]  to  analyze 
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(a)  A  simple  periodic  structure. 


(b)  A  compound  periodic  structure. 


Figure  1.1.  Examples  of  periodic  structures. 
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forced  response  of  such  periodic  beams.  Natural  frequencies 
and  principal  modes  were  also  found  for  periodic  beams  with 
certain  extreme  ends  by  Denke  et  al[5].  The  transfer  matrix 
method  was  shown  to  be  suitable  for  the  vibration  analysis 
of  such  multi-span  beams,  as  the  extreme  boundary  conditions 
can  be  easily  taken  into  account  with  the  matrix  equation. 

In  addition,  by  working  with  a  smaller  sized  matrix 
representing  a  single  periodic  beam  element  the 
computational  effort  becomes  virtually  independent  of  the 
number  of  spans.  For  years  the  transfer  matrix  method  has 
been  the  major  tool  for  vibration  analysis  of  periodic 
beams.  Tables  of  analytical  transfer  matrices  for  commonly 
appearing  periodic  beam  elements  are  listed  in  [18]. 

C.  Propagation  Constant 

When  a  periodic  structure  contains  many  supports,  so 
many  modes  must  be  included  in  the  analysis  that  the 
classical  normal  mode  or  transfer  matrix  methods  become  very 
cumbersome.  Under  this  circumstance,  it  is  convenient  to 
regard  the  structure  as  being  infinite  in  extent.  Discrete 
principal  modes  do  not  exist  for  infinite  systems,  instead 
the  differential  equations  governing  linear  vibration  of  the 
system  yield  free  flexural  waves  which  attenuate  or 
propagate  in  the  absence  of  damping  and  external  forces. 

With  some  damping  present,  the  forced  response  of  a  finite 
structure  even  with  a  few  spans  can  be  approximately 
obtained  by  assuming  the  number  of  spans  being  infinite  and 
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by  using  the  wave  approach.  This  method,  still  giving 
acceptable  results,  would  require  far  less  computation. 
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When  a  periodic  structure  is  vibrating  in  one  of  its 
free  waves,  there  exists  a  constant  ratio  between  the 
amplitude  of  motions  in  the  adjacent  periodic  elements  [3]. 
The  natural  log  of  this  ratio  of  the  amplitudes  was  first 
called  the  "propagation  constant"  by  Heckl[6].  The  flexural 
waves  are  then  expressed  as  continuous  functions  of 
directional  propagation  constants. 

In  the  works  of  Mead[ 1 1 , 1 2 , 1 3 ]  and  Sen  Gupta[19,20]  the 
propagation  of  the  flexural  waves  in  uniform  beams  were 
extensively  investigated  by  the  use  of  the  propagation 
constant.  The  propagation  constant  in  their  works  was 
obtained  from  the  receptances  of  beam  elements  [2]  and  it 
was  shown  that  the  forced  response  of  infinite  beams  are 
easily  obtained  by  this  propagation  constant  method.  Even 
in  the  case  of  beams  of  finite  number  of  spans  they  also 
managed  to  express  some  boundary  conditions  in  terms  of 
propagation  constants.  The  mono-coupled  system  of  a  uniform 
beam  on  knife-edge  supports  was  reconsidered  by  Sen 
Gupta[20]  and  it  was  shown  that  the  natural  frequencies  are 
graphically  obtainable  from  the  curve  of  propagation 
constants.  However,  in  the  case  of  multi-coupled  systems 
the  representation  of  the  extreme  boundary  condition  in 
terms  of  propagation  constant  causes  considerable  difficulty 
[13]  and  hence  the  propagation  constant  method  is  hardly 
used  in  finding  the  natural  frequencies  of  finite  span  beams 
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on  other  than  knife-edge  supports. 

In  the  present  study  the  propagation  constant  is 
recognized  as  the  natural  log  of  an  eigenvalue  of  the 
transfer  matrix.  By  doing  so,  the  useful  features  of  the 
propagation  constant  are  still  retained  thus  expanding  the 
limit  of  the  transfer  matrix  method. 

D.  Finite  Element  Method 

The  transfer  matrix  and  the  propagation  constant 
methods  have  been  based  on  the  analytical  expression  of  the 
normal  mode  solutions  to  the  differential  equations  of  beam 
vibration.  Since  these  analytical  solutions  are  restricted 
to  only  simple  beam  shapes  so  are  the  transfer  matrix  and 
the  propagation  constant  methods. 

Often  periodic  structures  in  engineering  are  so  complex 
that  a  uniform  beam  model  no  longer  represents  the  practical 
structure  satisfactorily  and  therefore  a  numerical  technique 
such  as  the  Finite  Element  Method  is  required  to  analyze  the 
dynamics  of  the  realistic  models  of  these  complex  periodic 
structures.  In  fact,  almost  any  shape  of  structure, 
periodic  or  non-periodic,  can  be  analyzed  by  this  numerical 
technique  with  the  aid  of  digital  computers.  However,  as 
the  number  of  spans  increases  the  accuracy  of  the  numerical 
results  deteriorates  while  the  computational  costs  grow 
enormously  showing  the  limit  of  the  application  of  straight 
numerical  methods  for  the  vibration  analysis  of  periodic 
structures.  Due  to  such  limitations  of  the  straight 


' 
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numerical  methods  the  transfer  matrix  method  has  gained 
renewed  interest  for  the  vibration  analysis  of  periodic 
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structures . 

In  the  works  of  Orris  and  Petyt  [17]  the  method  of 
obtaining  the  propagation  constant  from  the  impedance  matrix 
formulated  by  Finite  Element  Method  was  discussed. 
Independently,  Meirovitch  and  Engels [ 1 4 , 1 5 ]  showed  that  a 
matrix  equation  can  be  formed  from  the  same  impedance  matrix 
using  the  so  called  Z-transform  method.  Again,  by  defining 
the  propagation  constant  as  the  logarithmn  of  the  eigenvalue 
of  the  transfer  matrix,  these  two  seemingly  different 
techniques  can  be  combined  with  added  convenience  resulting 
from  the  properties  of  the  impedance  matrix. 

The  transfer  matrix  method  employed  in  the  present 
study  is  different  from  the  classical  one  in  that  the 
expression  is  derived  from  the  impedance  matrix  resulting 
from  the  Finite  Element  Method  not  from  the  normal  mode 
solution  to  differential  equations.  Thus  the  transfer 
matrix  for  a  substructure  is  obtainable  no  matter  how 
complex  it  may  be.  The  analytical  expression  for  the 
transfer  matrix  is  also  possible  as  a  special  case  if  the 
impedance  matrix  is  expressed  in  terms  of  the  normal  mode 
solutions  of  the  beam  differential  equations.  The 
analytical  expression  for  the  propagation  constant  is  then 
obtained  from  the  eigenvalues  of  the  analytically  derived 
transfer  matrix. 
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E.  Timoshenko  Beam 

All  the  analytical  studies  of  the  flexural  motions  of 
elastic  beams  resting  on  multiple  supports  have  been  based 
on  the  Bernoulli-Euler  beam  theory  neglecting  the  effect  of 
rotary  inertia  and  shear  deformations.  Corrections  to  the 
Bernoulli-Euler  beam  theory  were  given  by  Timoshenko [ 2 1 ]  in 
which  the  governing  differential  equation  includes  the 
effect  of  rotary  inertia  and  shear  deformation.  Huang[7] 
solved  this  differential  equation  by  forming  the  normal  mode 
equations  for  six  common  types  of  single  span  beams  and 
discussed  the  effect  of  cross-sectional  dimensions  on  the 
natural  frequencies  of  a  single  span  beam.  However,  little 
work  has  been  done  for  the  Timoshenko  beam  with  intermediate 
supports.  This  prompted  the  choice  of  a  multi-span 
Timoshenko  beam  as  the  numerical  example  for  the  present 
study . 

Quite  recently  several  Finite  Element  models  have  been 
proposed  for  the  same  Timoshenko  beam  as  considered  by 
Huang.  Kapur' s [ 8 ]  and  Craggs'[4]  elements  are  two  of  these 
that  have  given  good  results  over  a  wide  range  of 
cross-sectional  parameters.  These  Finite  Element  analyses 
as  well  as  Huang’s  analytical  study  of  the  beam  vibration 
show  that  when  solutions  are  based  on  the  Bernoulli-Euler 
beam  theory  there  are  considerable  differences  except  in  the 
low  frequency  range  of  a  long  slender  beam  where  the 
cross-sectional  dimension  is  negligible  with  respect  to  the 


span  length. 


. 


In  this  study  the  results  are  obtained  analytically  as 


well  as  numerically  for  the  Timoshenko  beam  on 
multi-supports.  These  results  are  then  compared  with  those 
based  on  Bernoull i-Euler  beam  theory  in  order  to  investigat 
the  effect  of  cross-sectional  dimensions  on  the  natural 
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II.  TRANSFER  MATRIX  AND  PROPAGATION  CONSTANT 


A.  Introduction 

When  a  periodic  structure  is  approximated  as  a  uniform 
beam  on  multiple  supports,  the  transfer  matrix  can  be  easily 
formulated  from  the  analytical  solutions  [9,10]  or  listed  in 
[18].  However,  many  engineering  structures  are  often  of  the 
form  which  may  not  be  approximated  as  a  uniform  beam,  and 
can  only  be  analysed  by  a  numerical  technique  such  as  the 
Finite  Element  Method.  This  requires  a  more  general  method 
of  deriving  the  transfer  matrix  for  any  type  of  periodic 
structure . 

In  the  following  section  the  transfer  matrix  which 
relates  the  state  vectors  at  two  successive  nodes  is  derived 
from  the  impedance  matrix  commonly  formulated  by  the  Finite 
Element  Method.  The  impedance  matrix  of  a  single  periodic 
unit  relates  the  intercoupling  generalized  forces  to  the 
generalized  displacements  at  both  ends  of  the  substructure. 
This  impedance  matrix,  expressed  either  analytically  or 
numerically,  displays  a  general  pattern  -  symmetry  is  one  of 
the  features.  Formulated  in  terms  of  submatrices  of  the 
impedance  matrix,  the  transfer  matrix  not  only  possesses  all 
the  properties  of  the  classical  transfer  matrix,  but  also 
contains  some  added  features  resulting  from  the  properties 
of  the  impedance  matrix.  For  this  reason  simple 
manipulations,  for  example  inversion,  of  the  transfer  matrix 
is  possible.  As  in  the  case  of  beams  supported  on 
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equi-spaced  elastic  springs,  the  change  of  configuration  of 
the  periodic  support  can  be  easily  done  by  simply  adding  or 
subtracting  the  appropriate  terms  to  the  diagonal  elements 
of  the  impedance  matrix. 

Once  the  transfer  matrix  is  formulated  and  its 
eigenvalues  are  found,  the  propagation  constants  are 
obtained  as  the  natural  logarithmns  of  these  eigenvalues. 

In  this  way  many  properties  of  the  propagation  constant  are 
simply  interpreted  by  the  algebraic  properties  resulting 
from  the  transfer  matrix  and  its  eigenvalues.  All  the 
considerations  of  the  behavior  of  transfer  matrices  and 
their  eigenvalues  can  be  interpreted  as  those  of  propagation 
constants  and  vice  versa.  Therefore,  only  a  brief 
discussion  of  the  propagation  constant  is  needed  in  this 
chapter . 

For  the  vibration  analysis  of  multi-unit  periodic 
structures,  either  free  or  forced,  often  the  multiplication 
of  the  element  transfer  matrix  derived  for  a  single 
substructure  is  required  in  order  to  find  the  total  transfer 
matrix  relating  the  state  vectors  from  one  end  to  another  of 
the  whole  structure.  This  direct  multiplication  is 
satisfactory  when  the  number  of  spans  is  small,  however,  for 
periodic  structures  composed  of  many  substructures  this 
direct  multiplication  would  be  very  costly  yielding  doubtful 
results.  In  this  circumstance,  either  the  Cayley-Hami lton 
theorem  or  the  similarity  transform  method  can  be  used  so 
that  the  accuracy  of  the  numerical  results  improve  while 
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reducing  the  computational  efforts. 

In  some  cases  only  the  propagation  constants  (or  the 
eigenvalues  of  the  transfer  matrix)  are  needed  instead  of 
the  transfer  matrix  itself.  Since  the  transfer  matrix  is 
expressed  in  terms  of  submatrices  of  the  impedance  matrix, 
the  eigenvalues  of  the  transfer  matrix  or  the  propagation 
constant  can  be  related  directly  to  the  impedance  matrix. 
In  this  way  much  of  the  computational  effort  is  reduced  by 
eliminating  the  formulation  procedure  of  the  transfer 
matrix . 

B.  Transfer  Matrix  from  Impedance  Matrix 

Consider  a  periodic  system  shown  in  Figure  2.1,  and 
suppose  that  the  impedance  matrix  [Z]  for  the  subsystem 


Figure  2.1.  Block  diagram  of  a  periodic  structure 

showing  the  generalized  displacements  and  forces. 
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shown  in  Figure  2.1  is  formed  as 

Cf }  =  [Z]{d]  (2.1) 

where  { d }  and  {f}  denote  the  generalized  displacement  and 
force  vectors,  respectively.  The  impedance  matrix,  [Z]  is 
of  the  familiar  form 


[Z]  =  [K]-w2[M] 


(2.2) 


when  formulated  by  the  Finite  Element  Method.  Here,  [K]  and 
[M]  denote  the  stiffness  and  mass  matrices,  respectively. 
This  impedance  matrix  as  well  as  the  displacement  and  force 
vectors  can  be  partitioned  as 


f 


L 


R 


~zLL 
_  Z  R  L 


(2.3) 


where  the  subscripts  L  and  R  denote  the  left  and  right  ends 
of  the  subsystem.  In  equation  (2.3),  it  is  assumed  that  the 
internal  nodal  degrees  of  freedom  are  already  removed,  if 
necessary,  by  the  procedure  described  in  [14], 

It  can  be  easily  shown  that  the  equation  (2.3)  can  be 
rearranged  in  the  form 
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But , 


{dR}  j  =  {d L }  j  +  ! 


(2.5) 


by  the  compatibility  relationship  of  displacements  and 


{  f  r  }  j  -  (filj+i 


(2.6) 


by  the  equilibrium  condition  when  there  are  no  external 
forces  applied. 

Equations  (2.5)  and  (2.6)  are  substituted  back  into 
equation  (2.4)  to  yield  a  matrix  difference  equation 
relating  generalized  displacements  and  forces  in  the 
successive  subsystems  j  and  j+1  as 


(2.7) 


[ T ]  denotes  the  transfer  matrix  by  definition  consisting  of 
submatrices  as 


~  Z  [  R  Z  L  L 


[T] 


(2.8) 


The  augmented  vector  of  generalized  displacements  and  forces 
is  called  the  state  vector  at  a  node. 
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C.  Properties  of  Transfer  Matrix 


Inversion 

It  can  be  noted  that  the  equation  (2.3)  may  be 
rearranged  in  the  different  order  from  equation  (2.4)  as 


Substitution  of  equations  (2.5)  and  (2.6)  into  equation 
(2.9)  yields  a  transfer  matrix  in  the  opposite  direction  to 
equation  (2.8),  i.e.,  the  inverse  of  the  transfer  matrix, 
[T]  as  defined  in  equation  (2.8). 


(2.10) 


Thus  the  inverse  of  transfer  matrix,  [T- 1 3  is  another 
transfer  matrix  which  consists  of  the  submatrices  as 


[T-  1  ] 


-Zr,1  Z 


RL  "  R  R 


-Zrl1 


_  —  Z  l  L  Zrl'Zrr  +  ZlR  Z  l  L  Zrl 


(2.11) 


Direct  multiplication  of  [ T ]  and  [T" 1 3  in  equations  (2.8) 
and  (2.11)  confirms  an  identity  matrix,  i.e., 


[T] [T_ ’ ]  =  [I ] . 


(2.12) 


' 
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This  derivation  of  [T" 1 ]  expressed  by  equation  (2.11)  does 
not  require  any  additional  calculations.  If  the  transpose 
of  [T" 1 ]  is  taken  as 


[T-t  ] 


Z  L  R  Z  RrZ  RL 


(2.13) 


and  noting  that  the  impedance  matrix  [Z]  is  symmetric,  i.e., 

Zll  =  ZuLf  Zrr  =  ZRr  and  Z  =  Z  l  r  •  (2.14) 

then  the  equation  (2.13)  can  be  simplified  as 

[T-T]  = 

Comparison  of  equation  (2.15)  with  equation  (2.8) 
immediately  suggests  an  algorithm  for  inversion  of  the 
transfer  matrix  which  is  expressed  by  equation  (2.8): 

(a)  Partition  the  transfer  matrix  into  four  submatrices, 

(b)  Switch  the  left  top  and  right  bottom  submatrices, 

(c)  Change  the  sign  of  right  top  and  left  bottom 
submatrices , 

(d)  Take  the  transpose  of  the  matrix  obtained  in  (c). 


ZrrZlr  Z  rl  ZprZlpJ  ZLL 


L  -Z 


1 

R 


Z  L  O  Z  L  L 


(2.15) 
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Eigenvalues  and  Eigenvectors 

Assuming  that  a  constant  ratio  (eigenvalue),  X  exists 


such  that 

{dp}  =  X{dL  } , 

(2.16) 

and 


{ £r }  =  -M f l  }, 

(2.17) 

and  by  subsituting  equations  (2.16)  and  (2.17)  for  {dR}  and 


{fR},  matrix  equation  (2.3)  can  be  rewritten  as 

{ f L }  =  [Zll  ] { d l }  +  X [ Z L r ] { dL } 

(2.18a) 

—  X { f  l }  =  [ZRL]{dL}  +  X[ZRR]{dL}. 

(2.18b) 

Equation  (2.18a)  is  multiplied  by  X,  then  added  to  equation 
(2.18b)  to  yield 


[  X  2  Z  L  R  +X  (  ZL  L  +Zr  R  )  +  Z  R  L  ]  {  d  L  }  =  0. 

(2.19) 

Equation  (2.19)  does  not  have  an  immediate  solution 

for  X , 

except  when  it  is  a  scalar  equation.  However,  together  with 
a  trivial  relationship 


{XdL}  =  X{dL] 


(2.20) 
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equation  (2.19)  may  be  augmented  as 


Equation  (2.21)  is  of  the  same  size  as  the  transfer  matrix, 
[T]  and  allows  the  use  of  ordinary  eigenvalue  subroutines  to 
find  X  values  and  the  corresponding  eigenvectors  { d L } .  The 
generalized  force  vector,  {fL}  is  then  found  by  the 
relationship 

{ f U  }  =  [ZLL]{dL]  +  [ Z  l  r ] {XdL } .  (2.22) 

The  augmented  vector  of  {dL}  and  {fL}  obtained  from 
equations  (2.21)  and  (2.22)  is  the  eigenvector  of  the 
transfer  matrix. 

It  can  be  easily  shown  that  the  transfer  matrix  and  its 
inverse  possess  an  identical  characteristic  equation  and 
thus  the  same  set  of  eigenvalues.  It  follows  then  that  the 
eigenvalues  of  the  transfer  matrix  must  be  pairwise 
reciprocals. 

Multiplication  of  transfer  matrices 

According  to  the  Cayley-Hami lton  theorem,  a  matrix 
itself  satisfies  its  own  characteristic  equation  and  thus 
the  power  of  a  matrix  of  size  2n  may  be  expressed  by  a 
linear  combination  of  any  2n  linearly  independent 
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polynomials  of  the  matrix, 

n 

[TN]  =  2  [c  jTj  +  djT" j ] .  (2.23) 

j  =  1 

This  equation  must  be  satisfied  by  the  inverse  of  [T]  as 
well 


n 

[T- N]  =  2  [c jT- j  +  djTj ] .  (2.24) 

J=1 


Addition  and  subtraction  of  equations  (2.23)  and  (2.24) 
yield 


n 

[Tn+t-N]  =  2  a  j [ T j  +T~  j ]  (2.25a) 

j=l 

n 

[ tn_t- n ]  =  2bj[Tj-T-j]  (2.25b) 

J=1 


where  n  denotes  half  the  size  of  the  transfer  matrix.  The 
constants  a j  and  b j  are  determined  by  substituting  the 
eigenvalues  X  for  [T]  in  equations  (2.25a)  and  (2.25b).  The 
A/-th  power  of  the  transfer  matrix,  [TN]  is  then  obtained 
f  rom 


n  n 

[T N]  =  2  a  j [T j  +T" j ]  +  2  bj[T'-T-j]  (2.26) 

i  =  1  j=i 

It  is  often  more  convenient  to  find  [Tn]  by  the 
so-called  similarity  transform  method.  Since  the  transfer 
matrix  has  2n  linearly  independent  eigenvectors  it  can  be 
diagonalized  as 


«> 
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[T]  =  [SHAHS-1]  (2.27) 

where  the  columns  of  the  matrix,  [S]  are  the  eigenvectors  of 
[T]  and  [A]  is  a  diagonal  matrix  with  the  corresponding 
eigenvalues  of  [ T ]  along  its  diagonal.  Then  the 
multiplication  of  [T]  can  be  done  as 

[TN]  =  [S][An][S'1 ].  (2.28) 


D.  Propagation  Constants  and  Transfer  Matrix 

Although  the  concept  of  the  propagation  constant  was 
originally  developed  in  wave  theory  [6],  independent  of  the 
transfer  matrix,  it  is  quite  clear  that  there  exists  a 
definite  relationship  between  them.  By  definition 

{v}j+1  =  [ T ] { v } j  :  transfer  matrix 

(2.29) 

{v}j  +  1  =  exp(At){v}j  :  propagation  constant 

where  {v}  denotes  the  state  vector  consisting  of  the 
generalized  displacements  and  forces  as  shown  in  equation 
(2.8),  and  u  is  called  the  propagation  constant. 

The  latter  can  be  subtracted  from  the  former  in 
equation  (2.29)  to  yield 


[T  -  exp(jx)  I  ]  { v }  j  =  0  . 


(2.30) 
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Thus  the  propagation  constant,  n  may  be  understood  as  the 
natural  logarithmns  of  the  eigenvalues  of  the  transfer 
matrix.  Hence  the  propagation  constant  in  this  study  does 
not  require  a  separate  formulation,  but  is  evaluated  from 
the  transfer  matrix  and  its  eigenvalues. 

The  propagation  constants  as  well  as  the  eigenvalues  of 
the  transfer  matrix  are  generally  complex  numbers.  Since 
the  eigenvalues  of  the  transfer  matrix  are  pairs  of 
reciprocals,  the  propagation  constants  are  pairs  of  opposite 
signs . 

For  the  transfer  matrix  of  size  2n ,  there  exist  n  pairs 
of  propagation  constants  and  hence  for  a  mono-coupled  system 
only  a  single  curve  is  needed  to  show  the  propagation 
constant . 

E.  Discussion 

The  transfer  matrix  obtained  from  the  impedance  matrix 
as  in  the  present  study  can  be  applied  to  any  complex 
periodic  structure  since  the  impedance  matrix  can  be  always 
formed  by  the  Finite  Element  Method.  The  Finite  Element 
Method  commonly  uses  polynomial  interpolation  functions,  for 
instance,  a  cubic  function  is  normally  assumed  for  the 
deflection  of  beams.  The  accuracy  of  the  Finite  Element 
Method  depends  mainly  on  how  closely  these  interpolation 
functions  approximate  the  exact  normal  modes.  If  these 
interpolation  functions  happen  to  be  identical  to  the  normal 
modes,  the  Finite  Element  Method  yields  the  exact  solution 
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within  the  round-off  error  due  to  the  digital  computation. 

In  the  case  of  a  uniform  beam  the  analytical  expression  for 
the  transfer  matrix  can  be  obtained  through  the  procedure  as 
described  in  the  section  II. B,  from  the  impedance  matrix 
analytically  formulated  using  the  two  hyperbolic  and 
trigonometric  normal  modes  as  the  interpolation  functions. 

The  cross-symmetric  property  of  the  transfer  matrix  has 
been  considered  quite  important  in  the  earlier  studies 
[9,10],  This  cross-symmetry  of  a  transfer  matrix  is 
expressed  in  the  indicial  form  as 

tjj=t2n+1-j,2n+1-i  (2.31) 


where  tjj  is  the  /- th  row  and  j- th  column  entry  of  the 
transfer  matrix  of  size  2n.  This  cross-symmetry  together 
with  other  properties  of  the  transfer  matrix  leads  to  an 
algorithm  for  transfer  matrix  inversion  as 


t ,  j  =  (-1)i  +  j  t  f  j  (2.32) 

where  tu  is  the  /-th  row  and  j-th  column  entry  of  [T' 1 ] . 

When  the  periodic  unit  is  symmetric  about  its  center  it 
is  possible  to  obtain  a  cross— symmetr ic  transfer  matrix  by  a 
suitable  choice  of  the  order  and  sign  convention  for  the 
components  of  the  state  vector.  However,  the  transfer 
matrices  are  not  generally  cross-symmetric,  and  even  so,  the 
proper  reordering  and  rearranging  of  a  transfer  matrix  into 
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a  cross-symmetric  form,  especially  with  the  digital  computer 
is  extremely  difficult.  In  the  present  study  dealing  with 
the  generally  shaped  periodic  structures,  the  use  of 
equation  (2.32)  can  not  be  made  and  thus  a  new  algorithm  of 
inversion  which  does  not  require  the  cross-symmetry  be 
developed  as  in  the  section  II. C. 

Certainly  the  eigenvalues  and  eigenvectors  can  be  found 
by  ordinary  subroutines  directly  from  the  transfer  matrix. 
When  its  application  is  restricted  to  the  periodic  beams 
(the  size  of  the  transfer  matrix  is  at  most  four)  it  matters 
little  how  the  eigenvalues  are  obtained.  However,  it  is 
possible  to  greatly  reduce  the  computational  effort  for 
finding  the  eigenvalues  and  eigenvectors  by  using  the 
submatrices  of  the  impedance  matrix  instead  of  the  transfer 
matrix  as  described  in  equations  (2.16)  through  (2.22). 

When  the  periodic  unit  is  symmetric  further  reduction  of 
computation  is  possible  by  halving  the  size  of  matrix 
involved  [14,15].  This  is  possible  because  of  the  special 
pattern  of  the  impedance  matrix  of  the  symmetric  subsystem, 
although  it  was  not  pointed  out  in  his  study.  By  halving 
the  size  of  matrix  for  the  eigenvalue  subroutine  the 
computational  cost  is  reduced  by  a  factor  of  approximately 
eight . 

By  expressing  the  transfer  matrix  in  terms  of  the 
submatrices  of  the  impedance  matrix  it  can  be  easily  shown 
that  the  transfer  matrix  and  its  inverse  possess  the  same 
characteristic  equation  and  thus  the  same  set  of 


. 
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eigenvalues.  Consequently,  the  determinant  of  a  transfer 
matrix  is  always  unity  and  the  eigenvalues  of  the  transfer 
matrix  are  therefore  pairwise  reciprocals.  This  leads  to 
the  conclusion  that  the  propagation  constants  which  are 
defined  as  the  logarithmns  of  the  eigenvalues  must  be  pairs 
of  opposite  signs.  Physically  this  should  be  true  since 
nature  would  not  recognize  the  negative  and  positive 
directions  of  the  structure. 

Without  excitation,  a  wave  signal  is  not  amplified  as 
it  propagates  along  the  structure.  This  implies  that  an 
eigenvalue  whose  magnitude  is  greater  than  unity  corresponds 
to  a  negative  going  wave  and  its  reciprocal  to  a  positive 
going  wave.  In  terms  of  propagation  constants,  the  one  with 
a  positive  real  part  corresponds  to  a  negative  going  wave 
while  the  one  with  a  negative  real  part  to  a  postive  going 
wave . 

Generally,  the  propagation  constant  as  well  as  an 
eigenvalue  of  the  transfer  matrix  is  a  complex  number.  The 
real  part  of  the  propagation  constant  represents  the 
amplitude  decay;  the  imaginary  part  the  phase  difference  at 
two  successive  nodes.  An  energy  carrying  free  flexural  wave 
is  possible  in  the  distinct  frequency  zones  where  the  real 
part  of  propagation  constant  vanishes.  For  a  finite 
structure  the  wave  motion  satisfies  the  extreme  end 
conditions  only  at  the  certain  groups  of  discrete  values  of 
frequencies  within  these  propagation  zones.  These  groups  of 
frequencies  constitute  the  natural  frequencies  of  the  whole 
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structure,  which  will  be  discussed  in  the  next  chapter. 

In  order  to  use  the  transfer  matrix  for  the  vibration 
analysis  of  the  A/-span  periodic  structure,  free  or  forced, 
the  first  step  required  may  be  the  multiplication  of  the 
element  transfer  matrix.  Three  different  techniques  are 
suggested,  namely;  direct  multiplication,  application  of 
Cay ley-Hami 1 ton  theorem  and  the  similarity  transform  method. 
Direct  multiplication  is  of  course  the  simplest  and  can  be 
satisfactorily  used  for  a  structure  with  a  few  spans  (three 
or  four  spans  at  the  most).  This  direct  multiplication  is 
the  only  available  technique  in  the  case  where  the  structure 
is  composed  of  different  substructures.  Considering  that 
the  matrix  multiplication  is  a  costly  process  losing  the 
accuracy  of  the  numerical  results,  this  direct 
multiplication  can  not  be  used  for  a  periodic  structure  of 
many  spans. 

The  use  of  the  Cayley-Hami 1 ton  theorem  also  contains 
the  direct  multiplication  2n  times  (n  times  for  [ T ]  and  n 
times  for  [T'1]  ).  Because  of  this  direct  multiplication 
involved,  this  method  is  not  suitable  for  multi-coupled 
systems  when  there  are  many  coupling  coordinates.  For  the 
sake  of  efficiency,  the  Cayley-Hami lton  theorem  for  the 
multiplication  of  the  transfer  matrix  is  not  directly 
applied  to  the  transfer  matrix  itself,  but  to  the  sum  and 
the  difference  of  the  transfer  matrix  and  its  inverse.  In 
this  way,  the  computation  deals  with  two  half-sized  matrices 
instead  of  one  double-sized  and  this  would  result  in 
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reduction  of  the  computational  effort.  Since  the  sum  and 
the  difference  of  the  transfer  matrix  and  its  inverse  are  of 
a  weakly  coupled  form,  with  half  the  components  being  equal 
to  or  close  to  zero,  the  accuracy  of  the  results  also 
improves . 

On  the  other  hand,  the  method  of  similarity  transform 
does  not  require  any  multiplication.  Instead,  it  requires 
that  the  eigenvectors  of  the  transfer  matrix  be  found  to 
form  the  matrix  [S]  and  its  inverse,  [S'1].  The 
eigenvectors  are  automatically  obtained  in  most  eigenvalue 
subroutines,  however,  the  inversion  of  [S]  is  again  a  costly 
process.  Fortunately,  many  vibration  problems  do  not  need 
[S'1]  to  be  found  as  it  will  be  shown  in  the  next  chapter. 

In  such  cases  the  method  of  similarity  transform  can  be  most 
efficiently  used. 


■: 


III.  HARMONIC  MOTIONS  OF  PERIODIC  STRUCTURES 


A.  Introduction 

In  the  previous  chapter  a  transfer  matrix  was 
formulated  for  a  typical  bay  of  a  periodic  structure.  This 
element  transfer  matrix  was  then  shown  to  be  multiplied  to 
represent  a  transfer  matrix  for  an  N~ bay  periodic  structure, 


i  .  e 
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(a ) .  In  this  case 

the 

t ransf e 

r  matrix  f 

or  the 

tw 

o  end 

bays  may  be  formulated  separately  and  the  transfer  matrix 
for  the  whole  system  can  be  obtained  as 


{v}N  =  [T]  [TN'2]  C  T ] ,  {v}0  (3.2) 


where 

[T]  , 

and 

[ T ] N  denote 

the 

first 

and 

last 

elements  of 

the 

many  homogeneous  boundary  beam 
be  considered  identical  to  the 


transfer  matrices  for  the 
structure.  Fortunately  in 
cases,  the  two  end  bays  can 
other  interior  bays.  These 
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are  the  cases  where  the  extreme  boundaries  offer  more 
constraints  than  the  intermediate  supports  -  see  Figure 
3.1(b).  The  boundary  condition  is  then  applied  to  the  state 
vectors  {v}0f  {v}N  in  equation  (3.1)  to  obtain  the  frequency 
equation  of  the  complete  system.  Therefore,  without  loss  of 
generality,  it  is  assumed  in  the  present  study  that  [TN] 
represents  the  transfer  matrix  for  the  whole  structure 
relating  the  state  vectors  {v}0  and  {v}N. 

Special  consideration  is  given  to  mono-coupled  systems 
where  the  periodic  units  are  inter-coupled  by  a  single 
displacement  or  a  force  (generalized).  A  spring-mass 
oscillating  system  and  a  uniform  beam  resting  on  knife-edge 
supports  are  the  examples  of  a  mono-coupled  system.  The 
transfer  matrix  for  a  mono-coupled  system  is  of  size  two, 
and  hence  it  has  a  pair  of  eigenvalues  and  consequently  a 
pair  of  propagation  constants.  This  propagation  constant 
can  be  represented  by  a  single  curve  and  it  is  possible  to 
use  a  graphical  technique  to  find  the  natural  frequencies  of 
the  system.  In  the  case  of  a  multi-coupled  system,  the  set 
of  frequency  equations  is  of  determinantal  form,  consisting 
of  implicit  functions  of  frequencies,  which  can  be  solved  by 
an  iterative  technique  or  a  numerical  scheme  which  does  not 
require  explicit  knowledge  of  frequencies. 

While  it  is  possible  to  consider  transient  or 
non-harmonic  excitation  problems,  only  the  steady-state 
harmonic  forced  response  problems  are  considered  in  this 
study.  In  such  cases  the  response  can  be  uniquely 
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(a)  End  transfer  matrices 
are  required. 


(b)  End  transfer  matrices 
are  not  required. 


Figure  3.1.  Uniform  beams  with  end  elements  different 
from  the  periodic  element. 
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determined  without  any  iterative  procedure. 

When  a  periodic  structure  contains  infinitely  many 
spans  the  frequency  equation  is  no  longer  obtained. 

However,  the  steady-state  solutions  to  the  harmonic 
excitation  are  still  available  from  relatively  simple 
calculations.  When  this  response  is  shown  in  the  form  of  a 
resonance  curve  the  frequency  at  which  the  resonance  occurs 
can  be  considered  as  the  natural  frequency  of  such  an 
infinite  system. 

B.  Natural  Frequencies 

Consider  a  periodic  structure  consisting  of  N  periodic 
substructures  with  2n  degrees  of  freedom  at  each  node. 

There  exist  2n  boundary  conditions  -  n  for  the  left  boundary 
and  n  for  the  right.  These  boundary  conditions  are 
substituted  into  equation  (3.1)  and  a  set  of  equations  for 
the  boundary  conditions  is  obtained.  If  the  boundary 
condition  is  such  that  .all  the  generalized  displacements  are 
zeros  in  the  state  vectors,  {v}0  and  { v } N  so  that 

[T,n2]  {  f  }  o  =  0  (3.3) 

where  [T^]  denotes  the  top-right  submatrix  of  [TN]. 

Equation  (3.3)  represents  a  set  of  n  homogeneous  algebraic 
equations,  which  has  nontrivial  solutions  if  and  only  if  the 
determinant  of  [T^]  vanishes 
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det  [T,n2  ]  =  0.  (3.4) 

The  frequency  equation  can  be  obtained  differently  if 
the  similarity  transform  method  is  used  for  obtaining  [TN]. 
From  equation  (2.28) 

{ v } n  =  [S][AN][S-’]  {v} o .  (3.5) 

Premultiplication  of  both  sides  by  [S-1]  yields 

[S-'Uv}H  =  [ AN ] [ S "  1  ]  { v  }  o  -  (3.6) 


By  setting 


[S]{C}n={v}n/  (3.7) 

and 

[S]  {C} 0  =  {v}0,  <3-8> 

equation  (3.5)  can  be  rewritten  as 

{C}N  =  [AN]  { C } o •  (3*9) 

In  these  equations  { C } 0  and  {C}N  denote  the  arbitrary 
coefficients  to  be  determined  by  the  boundary  conditions  at 
the  left  and  right  ends  of  the  structure.  Applying  the 

conditions  to  equations  (3.7)  and  (3.8)  a 


extreme  boundary 
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set  of  2n  homogeneous  frequency  equations  are  obtained.  If 
the  both  ends  are  fixed,  i.e.,  all  the  generalized 
displacements  are  zeros,  then 


[s]  !  {  C  }  0  =  0 
[SAN]  ,  (do  =  o 


(3.10) 


where  the  subscript  1  in  equation  (3.10)  denotes  the  upper 
half  of  the  matrix  and  the  total  number  of  equations  in 
(3.10)  is  2n ,  Boundary  conditions  other  than  fixed  ends  can 
be  handled  similarly.  These  frequency  equations  (3.4)  or 
(3.10)  are  highly  transcendental  and  not  to  be  solved 
simply,  because  the  components  of  the  matrices  [Tn]  and  [S] 
are  implicit  functions  of  frequencies.  Thus  a  numerical 
algorithm  which  does  not  require  explicit  knowledge  of  the 
characteristic  determinant  is  needed  to  solve  the  frequency 
equation.  Meirovitch  and  Engels[15]  suggested  an  algorithm 
which  is  referred  to  as  Muller's  method.  For  the  detailed 
discussion  for  the  Muller’s  method  one  should  refer  to  [15]. 

C.  Steady-state  Harmonic  Response 

Consider  again  an  A/-bay  structure  and  suppose  that  a 
harmonic  excitation  is  introduced  at  the  node  j.  The  state 
vectors  at  the  ends,  {v}0,  and  {v}N  are  now  related  as 

{ V }  n  =  [T"l(v)o  +  [TN-JHw}j  (3.11) 
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where  {w} j  is  the  input  vector  at  j.  Equation  (3.11) 
contains  total  of  4n  unknown  quantities  in  2n  equations. 

Upon  imposing  2n  boundary  conditions  on  {v}0  and  {v}N,  the 
number  of  total  unknowns  is  reduced  to  2n ,  giving  a  set  of 
2n  equations  with  2n  unknowns.  This  equation  is 
inhomogeneous  due  to  the  input  vector  and  thus  the  state 
vectors  {v}0  or  {v}N  can  be  uniquely  determined.  The 
response  at  the  intermediate  node,  /  then  can  be  obtained  by 


{ v } j  =  [ T ' ] { v } o :  if  i^ j , 

{v} ;  =  [ T ' ] { v } o  +  [Ti_j]{w}j:  if  i>j. 


(3.12) 


In  the  case  where  the  input  vector  is  introduced  in 
between  the  nodes  a  partial  transfer  matrix  is  needed,  while 
the  rest  of  the  analysis  is  the  same.  Suppose  an  excitation 
is  introduced  at  a  point  x  in  between  the  nodes  j- 1  and  j, 
then  the  equation  (3.11)  can  be  modified  as 

{v} N  =  [ T  N ] { v } o  +  [TN-j][Tj.x]{w}x  (3.13) 

where  [Tj_x]  denotes  a  partial  transfer  matrix  transferring 
the  input  vector  from  the  point  x  to  the  node  j.  When  there 
are  multiple  input  vectors  the  term  due  to  other  input 
vectors  can  be  simply  added  to  the  RHS  of  equation  (3.11). 

The  problem  of  the  periodic  structure  consisting  of 
infinitely  many  spans  is  not  difficult  to  solve  but  needs 
separate  consideration.  Figure  3.2(a)  shows  an  example  of  a 
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semi- inf  ini te  system.  The  boundary  condition  at  the  right 
end  is  that  there  is  no  wave  returning  from  it  so  that  half 
the  coefficients  {C}0  in  equations  (3.7)  and  (3.8)  which 
correspond  to  negative  going  waves  are  zeros.  This  gives  n 
boundary  conditions  together  with  another  n  boundary 
conditions  at  the  other  end  provide  2n  required  equations. 


W 

(a)  Semi-infinite  system. 


W 

(b)  Infinite-infinite  system. 


Figure  3.2.  Harmonically  excited  infinite  systems. 


Consider  an  example  of  a  semi-infinite  system  with  the 
excitation  at  the  free  end.  From  equation  (3.6) 

{v} o+{w} o  =  [ S ] { C } o  (3.14) 


where  half  the  elements  of  {C}0  which  correspond  to  negative 
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going  waves  are  zeros  and  the  bottom  half  of  (v}0  are  also 
zeros  so  that  n  equations  from  the  bottom  are  solved  first 
for  the  n  nonzero  coefficients  then  the  response  is  given  by 
the  upper  half  of  the  equations.  In  the  case  of  the 
infinite-infinite  system  with  the  exciting  force  introduced 
at  the  middle  as  shown  in  Figure  3.2(b),  the  additional  n 
equations  are  the  conditions  at  the  excitation  that  the  n 
generalized  displacements  are  continuous. 

D.  Mono-coupled  System 

Figure  3.3  shows  two  examples  of  mono-coupled  system. 
The  2x2  transfer  matrix  for  a  mono-coupled  system  is  of  the 
form 


[T] 


til  t  i  2  1 

-  t  2  1  t  2  2  J  • 


(3.15) 


Since  the  determinantal  value  of  a  transfer 
unity  the  inverse  of  [T]  in  equation  (3.15) 


matrix  is  always 
is  simply 


[T-  1  ] 


(3.16) 


Let  the  eigenvalues  of  [T]  be  X  and  1/X  and  the  propagation 
constant  be  ±m,  then  the  total  transfer  matrix,  [T  ]  is 
obtained  by  the  equations  (2.25)  and  (2.26)  as 


<  ' 
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[Tn] 


coshNM+ 


sinhNM 

sinh*/ 


sinhN^t  . 
sinh*t  tzi 


lt  i  i  ~t  2  2 )  s  inhN/i 
2  sinhji  tl2 

coshN<J-^^-it.i'-t2Z) 
sinhAi  2 


(3.17) 


where  it  is  noted  that 


coshu  =  0 . 5 ( X+ 1 /X ) 

-  0  •  5  (  t  i  1  ^"t  2  2  )  • 


(3.18) 


(a)  A  spring-mass  oscillating  system, 
m  m  2m 

*-Mh+  •— A/W — 


(b)  A  uniform  beam  on  knife-edge  supports. 


r 

V£> 

v£) 

A  A 

A 

A 

A 

A 

Figure  3.3.  Examples  of  mono-coupled  systems  showing 
two  possible  element  divisions. 
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In  the  special  case  where  the  periodic  element  is  symmetric 
about  its  center,  t i , =t 2 2 =cosh^ ,  so  that  equation  (3.17)  can 
be  simplified  as 


[TN] 


coshNju 

sinhNg 
-  sinh/i  2  1 


sinhNju . 
sinhM  12 

coshN/i 


(3.19) 


Three  different  boundary  conditions  may  be  considered; 

1 )  Both  ends  free 

sinhNM/sinh*i  =  0,  or  *z  =  j7ri/N  (  j=  1 , 2  ,  .  .  .  ,  N-  1  ) 
or  1 2 i  =  0 

2)  Both  ends  fixed 

sinhNM/sinhM  =  0  ,  or  a/  =  j 7r i /N  (  j=  1 , 2  ,  .  .  .  , N-  1  ) 
or  t12  =  0 

3)  One  end  free,  one  end  fixed 

1 .  Symmetric  element 

coshN^  =  0,  or  n=(  j  +  0.5)7ri/N  (  j  =  0  ,  1  ,  .  .  ,  N-  1  ) 

2.  Unsymmetric  element 

coshN^iO.  5sinhNMt  ii“t22)/sinh/i  =0. 

For  the  spring-mass  oscillating  system  shown  in  Figure 
3.3(a),  the  transfer  matrix  is  obtained  as  -  see  Appendix 


(3.20) 


(3.21 ) 


(3.22) 


(3.23) 


. 


38 


[T] 


1-b2  -1  ' 

- b 2 ( 2-b 2 )  1-b2J, 


(3.24) 


where  b2 -mo;2 /k .  The  propagation  constant  is  obtained  from 
this  transfer  matrix  as 


M  =  cosh" 1 ( 1-b2 )  (3.25) 

which  is  shown  in  Figure  3.4.  Above  the  frequency,  b=/2  the 
wave  attenuates  with  increasing  amplitude  decay  as  the 
frequency  increases  and  the  wave  propagation  is  only 
possible  in  the  frequency  zone,  b<i/2 ,  where  the  real  part  of 
the  propagation  constant  vanishes.  Within  the  propagation 
zone  at  the  points  numbered  1,2,...,  where  u=jni/U  (N=4  in 
Figure  3.4),  the  natural  frequencies  of  either  the  free-free 
or  the  fixed-fixed  systems  are  read,  i.e., 

b  =  / 1  -cos  (  j  7t/N  )  (j  =  1  ,2,  .  .  ,N— 1  ) 
and  additional  roots  exist 

b  =  0  and  /2 

in  the  case  of  free-free  boundary. 

The  natural  frequencies  of  the  system  with  the 
free-fixed  boundaries  fall  at  the  middle  of  these  points 
corresponding  to  the  natural  frequencies  of  the  free-free  or 
fixed-fixed  systems  as 


(3.26) 


(3.27) 


PROPAGATION  CONSTANT 
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Figure  3.4.  Propagation  constant  for  a  spring-mass  system. 
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b  =  i/l-cos(  (  j  +  0.5ir)/N)  (  j  =  0,  1  ,2,  .  .  ,N-1 )  (3.28) 

only  when  the  periodic  element  is  symmetric.  For  the  four 
span  spring-mass  systems  of  symmetric  elements  shown  in 
Figure  3.5  the  natural  frequencies  are  obtained  by  this 
method  and  they  are  given  in  Table  3.1.  The  mode  shapes 
corresponding  to  these  natural  frequencies  of  the  free-free 
and  fixed-fixed  cases  are  shown  in  Figure  3.6.  The  axial 
motion  of  the  mass  is  shown  vertically  in  this  figure. 

The  same  spring-mass  system  can  be  divided  into 
unsymmetric  periodic  elements  as  shown  in  Figure  3.3(a). 

The  transfer  matrix  in  this  case  becomes  -  see  Appendix  A. 2 


[T] 


"  1  -2b 2  -1  ‘ 

.  2b2  1  J. 


(3.29) 


The  propagation  constant  is  identically  obtained  as  in 
equation  (3.25)  and  Figure  3.4.  Figure  3.7  shows  the  four 
span  spring-mass  systems  of  unsymmetric  elements.  Note  that 
the  mass  at  the  ends  becomes  2m  in  this  case.  The  natural 
frequencies  of  the  symmetric  boundary  systems,  (a)  free-free 
and  (d)  fixed-fixed,  may  be  determined  graphically  from  the 
propagation  constant  curve  as  in  the  case  of  symmetric 
element  systems.  Whereas,  those  of  the  unsymmetric  boundary 
systems,  (b)  free-fixed  and  (c)  fixed-free,  can  be  obtained 
by  solving  the  frequency  equation  (3.23).  These  results  are 

given  in  Table  3.2. 


Table  3.1.  Natural  frequencies  of  four  span  spring-mass 

systems  of  symmetric  elements.  (b=J^cj) 


No.  MODE 

FREE-FREE 

FREE-FIXED 

FIXED-FIXED 

1 

o 

• 

o 

0.276 

0.541 

2 

0.541 

0.786 

1.0 

3 

1.0 

1.176 

1.307 

4 

1.307 

1.387 

5 

1.414 

m  k  m 


(a)  Free-free  boundary. 

m  2m  2m  2m  m 


(b)  Free-fixed  boundary. 

m  2m  2m  2m 

^ — • 


(c)  Fixed-fixed  boundary. 

2m  2m  2m 


Figure  3.5.  Free  vibrations  of  a  four  span  spring-mass 
system  of  symmetric  elements. 
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rain 


(b)  Fixed-fixed  system. 


(a)  Free-free  system. 


Figure  3.6.  Mode  shapes  of  the  four  span  spring-mass 


system . 
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Table  3.2.  Natural  frequencies  of  four  span  spring-mass 

systems  of  unsymmetric  elements.  (b=J^u)) 


MODE 

FREE-FREE 

FREE-FIXED 

FIXED-FREE 

FIXED-FIXED 

1 

o 

• 

o 

0.246 

0.315 

0.541 

2 

0.541 

0.707 

0.882 

1  .0 

3 

1  .0 

1  .083 

1  .274 

1  .307 

4 

1  .  307 

1  .  329 

5 

1.414 

2m  k 


(a)  Free-free  boundary. 

2m  2m  2m 


2m 


III* 
•  •  •  • 


(b)  Free-fixed  boundary. 


(c)  Fixed-free  boundary 


2m  2m  2m 


(d)  Fixed-fixed  boundary. 

2m  2m  2m 

— AW— 


Figure  3.7. 


Free  vibrations  o£  a  four  span  spring-mass 


system  of  unsymmetric  elements. 
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Next,  consider  the  spring-mass  systems  with  a  unit 
harmonic  force  as  shown  in  Figure  3.8.  For  the  first 
example  of  the  four  span  free-free  system  the  equation 
(3.11)  can  be  written  as 


(3.30) 


From  this  equation  the  displacement  of  the  first  mass,  X0 
can  be  obtained  as 


X0  -  ~  1 2  2 /t  2  1 

=  cosh4M • s i nh^/s inh4y ( 2b 2 -b M . 


(3.31  ) 


The  displacement  of  the  first  mass  in  the  next  example  of 
free-fixed  system  in  Figure  3.8  can  be  similarly  obtained  as 


Xo  -  ~  t i 2 /t 1 1 

=  sinh4/i/si.nh/u»cosh4/i. 


(3.32) 


The  displacement  of  the  center  mass  in  the  fixed-fixed  four 
span  spring-mass  system  can  be  also  obtained.  Since 


F0  can  be  obtained  as 


*; 
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m  k  m 
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#  -AW  #  AW  #  -VvW  -0  v\W  # 


Figure  3.8.  Forced  responses  of  symmetric  spring-mass 
systems . 
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F0  =  -sinh2;i/sinh4M . 

The  displacement  of  the  center  mass,  X2  is  then 

X2  =  t12F0 

=  sinh22/Lt/sinh|i*  sinh4M 


(3.34) 


(3.35) 


where  the  propagation  constant,  n  in  these  equations  can  be 
evaluated  by  the  equation  (3.25).  For  the  last  example  of 
infinite  system  in  Figure  3.8  the  response  can  be  obtained 
by  the  simpler  calculation.  Since  the  eigenvector  of  the 
transfer  matrix  in  equation  (3.24)  is  of  the  form 

X/F  =  1 / ±/b 4 -2b 2  (3.36) 


the  ratio  of  displacement  to  force  is  the  same  for  all  the 
state  vectors  as  in  this  equation.  The  harmonic  responses 
for  the  four  span  spring-mass  systems  in  Figure  3.8  obtained 
as  in  the  equations  (  3 . 3  1  ) - ( 3 . 36 )  are  shown  in  Figure 
3 . 9  ( a  )  -  ( d )  . 

E.  Discussion 

The  transfer  matrix  obtained  in  the  previous  chapter 
for  a  single  periodic  unit  of  the  periodic  structure  and  the 
propagation  constant  obtained  from  it  is  used  to  solve  the 
free  or  forced  vibration  problems  of  the  N  bay  structure  in 
this  chapter.  For  this  the  first  step  required  is  to 
multiply  the  element  transfer  matrix  in  order  to  obtain  the 
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(a)  Free-free  four-span  system. 


)  Free-fixed  four-span  system. 


Figure  3.9.  Resonance  curves  of  four  span  spring  mass 
systems . 
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(c)  Fixed-fixed  four-span  system. 


1.0  2.0 
FREQUENCY  b 


FREQUENCY  b 


Figure  3.9.  continued... 
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total  transfer  matrix  relating  the  state  vectors  at  the 
nodes  of  several  spans  apart,  eventually  from  one  end  to 
another.  It  is  shown  that  the  use  of  the  Cay ley-Hami 1 ton 
theorem  can  be  made  for  the  mono-coupled  system  quite 
effectively . 

Many  subroutines  for  the  eigenvalues  do  not  require 
much  additional  computation  to  find  the  eigenvectors.  These 
eigenvectors  are  simply  collected  to  form  the  matrix  [S]  and 
as  in  the  case  of  equation  (3.10)  the  inversion  of  this 
matrix  is  not  always  required.  That  is,  no  costly 
computation  with  the  matrices  is  involved  in  such  cases. 

This  also  increases  the  accuracy  of  the  numerical  results 
because  most  round-off  errors  occur  during  the  matrix 
multiplication  and  inversion  process.  The  elaborate 
numerical  scheme  of  avoiding  these  calculation  errors  as  in 
the  reference  [10]  is  not  necessary  in  the  present  study. 

The  natural  frequencies  of  the  periodic  structure  are 
found  where  the  determinant  of  a  submatrix  vanishes  not  by 
the  responses  becoming  infinite  as  in  [10].  Indeed  it  would 
be  a  difficult  task  to  obtain  accurate  results  with 
calculations  involving  infinity.  The  determinantal 
frequency  equations  (3.4)  and  (3.10)  can  be  solved  by  any 
numerical  techniques  available  for  such  equations.  However, 
it  was  found  that  the  determinant  behaves  very  smoothly 
around  the  roots  so  that  even  the  linear  interpolation  is 
quite  satisfactory.  On  the  other  hand  the  problem  of  the 
steady-state  response  to  harmonic  excitation  does  not 
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require  any  iterative  numerical  scheme,  instead  the 
solutions  are  obtained  by  a  linear  equation  solver  such  as 
Gaussian  elimination. 

The  forced  response  section  is  only  briefly  discussed 
in  this  chapter  because  there  are  numerous  types  of 
problems,  many  of  which  are  out  of  the  scope  of  the  present 
study.  It  is  believed  that  most  engineering  problems  of 
forced  responses  can  be  satisfactorily  solved  by  using  the 
superposition  of  harmonic  responses  or  by  the  classical 
normal  mode  method  using  the  results  obtained  by  the  present 
method . 

In  order  to  treat  the  infinite  structures  rigorously 
the  velocities  of  the  waves,  group  and  phase  velocities 
should  be  considered.  However,  the  infinite  structure  in 
this  study  is  not  really  infinite,  instead  it  approximates  a 
structure  with  many  spans  (more  than  five)  with  some  damping 
in  order  to  reduce  the  computational  effort  still  obtaining 
reasonably  good  results. 

An  entire  section  is  devoted  to  the  mono-coupled  system 
because  of  easy  calculation  and  easy  illustration  of  the 
method  and  the  behaviour  of  the  solutions.  The  propagation 
constant  of  the  mono-coupled  system  can  be  either  pure 
imaginary  or  real  with  the  imaginary  part  being  0  or  ±7 r,  in 
the  attenuation  zone.  For  a  spring-mass  system  only  the  low 
frequency  (  b  <  /2  )  signals  can  propagate  along  the  system 
so  that  the  system  acts  like  a  low  pass  filter.  On  the 
other  hand  uniform  beam  on  knife-edge  supports  acts  somewhat 
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like  an  alternating  frequency  band  pass  filter. 

For  an  N  bay  discrete  system  with  symmetric  units, 
there  exist  either  N  or  N~ 1  natural  frequencies  in  this 
propagation  zone,  with  the  interval  of  7r/A/  of  the  imaginary 
part  of  the  propagation  constant.  A  system  with  the  mixed 
boundary  has  N  masses  and  N  natural  frequencies  all  inside 
the  propagation  zone.  The  free-free  system  has  N+ 1  masses 
and  the  same  number  of  natural  frequencies,  of  which  N~ 1 
frequencies  lie  inside  the  propagation  zone  and  two  are  the 
bounding  frequencies.  The  lower  bounding  frequency  is  zero 
corresponding  to  rigid  body  motion  leaving  N  natural 
frequencies.  The  fixed-fixed  system  has  N~ 1  masses  and  N~ 1 
natural  frequencies  all  inside  the  propagation  zone. 

An  N  bay  continuous  system  has  N  natural  frequencies  in 
each  propagation  zone.  All  the  N  natural  frequencies  lie 
inside  the  propagation  zone  in  the  case  of  mixed  boundary 
system;  while  the  systems  with  free-free  or  fixed-fixed 
boundaries  have  N~  1  natural  frequencies  inside  the 
propagation  zone  together  with  one  of  the  bounding 
f  requenc ies . 

This  graphical  technique  of  finding  the  natural 
frequencies  was  first  noted  by  Sen  Gupta  [20]  by  considering 
the  flexural  waves  along  the  uniform  Bernoull i -Euler  beam 
resting  on  the  knife-edge  supports,  which  is  an  example  of 
the  mono-coupled,  continuous  system  as  shown  in  Figure 
3.3(b).  Figure  3.10  shows  the  propagation  constant  for  such 
beams  with  no  rotational  springs  attached  and  his  method  of 
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determining  the  natural  frequencies  of  the  four  span 
hinged-hinged  beams.  Unlike  the  spring-mass  system,  the 
attenuation  and  propagation  zones  alternate.  There  is  a 
group  of  N  natural  frequencies  in  each  propagation  zone,  and 
there  are  infinitely  many  such  groups  resulting  in 
infinitely  many  natural  frequencies.  This  method  of  finding 
the  natural  frequencies  of  the  uniform  Bernoull i-Euler  beams 
on  knife-edge  supports  is  identical  to  that  employed  for  the 
mono-coupled  systems  in  the  previous  section.  In  his  study, 
however,  the  limitations  of  this  graphical  technique  are  not 
clear.  In  the  present  study  the  mono-coupled  system  is  just 
a  special  case  where  the  present  method  of  analysis  applies 
and  in  the  procedure  of  the  analysis  the  graphical  technique 
appears.  It  is  now  clear  that  this  graphical  technique  can 
be  applied  to  any  mono-coupled  system,  whether  it  is  a 
uniform  beam  (Timoshenko  or  Bernoul 1 i -Euler  beam)  on 
knife-edge  supports,  spring-mass  oscillation  system  or  even 
some  electrical  circuits.  When  the  periodic  element  is  not 
symmetric  about  its  center,  however,  the  natural  frequencies 
of  the  system  with  unsymmetric  boundaries,  i.e.,  one  end 
free  and  the  other  fixed,  can  not  be  obtained  graphically, 
instead  the  equation  (3.23)  is  to  be  used. 

It  can  be  noted  that  the  symmetric  and  the 
anti-symmetric  modes  alternate  in  the  case  of  symmetric 
systems  in  Figure  3.6.  At  each  of  the  natural  frequencies 
of  the  free-free  system  including  the  first  one  which 
represents  the  rigid  body  motion,  the  end  masses  must 
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Figure  3.10.  Propagation  constant  for  a  Bernoull i-Euler 
beam  resting  on  knife-edge  supports. 
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oscillate  so  that  if  the  excitation  is  introduced  at  the  end 
the  resonance  must  occur.  However,  if  the  center  mass  of 
the  f ixed-f ixed  system  is  excited,  the  resonance  can  occur 
only  at  the  first  and  third  natural  frequencies 
corresponding  to  the  symmetric  modes.  This  may  be  observed 
in  Figures  3.9  (a)  and  (c). 

Harmonically  excited  steady-state  responses  are 
obtained  in  the  form  of  receptance  (ratio  of  displacement  to 
force  [2])  for  the  four  span  spring-mass  systems  with  four 
different  boundary  conditions  and  the  results  are  shown  in 
the  resonance  curves,  Figure  3.9(a)-(d).  Unlike  the  single 
mass  system  the  response  can  become  zero,  that  is,  a  force 
can  be  exerted  into  the  system  with  no  motion  at  the  point 
of  excitation;  the  system  is  acting  like  a  vibration 
absorber.  Since  the  curves  show  only  the  absolute  value  of 
the  displacement  to  force  ratio  the  phase  difference  of 
input  and  output  is  not  seen  from  these  figures  instead  it 
can  be  noted  from  the  equations  ( 3 . 3 1 ) - ( 3 . 36 ) .  The 
receptance  in  equation  (3.36)  is  pure  imaginary  when  b<]/2  so 
that  there  exists  0.5tt  phase  difference  in  the  infinite 
system.  The  free-free  system  is  capable  of  rigid  body 
motion  so  that  the  receptance  is  infinite  at  zero  frequency 
(static  load),  but  as  soon  as  the  exciting  force  becomes 
harmonic  with  low  frequency  the  displacement  becomes 
harmonic  with  i r  out  of  phase  of  the  force. 
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IV.  NUMERICAL  EXAMPLE 


UNIFORM  TIMOSHENKO  BEAM 


A.  Introduction 

All  the  formulation  in  the  previous  chapter  is  not 
restricted  to  simple  types  of  uniform  beams.  However,  the 
numerical  example  taken  in  this  chapter  is  a  uniform 
Timoshenko  beam  in  which  the  effect  of  rotary  inertia  and 
shear  deformation  is  fully  taken  into  account. 

In  this  chapter  the  impedance  matrix  for  a  beam  element 
is  first  obtained.  This  impedance  matrix  can  be  obtained  by 
the  Finite  Element  formulation  using  Huang’s  normal  modes  as 
the  interpolation  functions.  However,  the  same  impedance 
matrix  is  more  readily  available  by  the  beam  flexural  theory 
and  the  unknown  coefficient  method.  The  analytical  transfer 
matrix  for  a  Timoshenko  beam  is  then  derived  from  this 
impedance  matrix  and  expressed  in  terms  of  Huang's  normal 
mode  solutions  and  his  notations. 

Natural  frequencies  of  the  uniform  Timoshenko  beam  on 
the  knife-edge  supports  are  found  by  a  graphically 
illustrated  technique  as  explained  in  section  III.D, 
although  the  actual  calculation  is  done  by  a  digital 
computer.  The  natural  frequencies  of  beams  on  spring 
supports  are  found  by  solving  the  determi nan tal  frequency 
equations.  The  configuration  of  the  intermediate  supports 
is  changed  by  varying  the  translational  spring  constant,  /<t 
and  the  rotational  spring,  kr.  For  example,  the  beam  resting 
on  the  knife-edge  supports  can  be  represented  either  by  a 
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beam  with  kt  becoming  infinite,  or  by  eliminating  the 
deflection  and  the  shear  force  at  the  nodes  since  only  the 
bending  slope  and  the  bending  moment  can  be  transmitted 
through  the  knife-edge  supports. 

The  rotational  spring  in  Figure  4.1  is  assumed  to  exert 
a  resisting  moment  proportional  to  the  bending  slope  instead 
of  the  total  slope.  Consequently,  the  boundary  condition  at 
the  clamped  end  of  the  Timoshenko  beam  is  implemented  by 
equating  the  bending  slope  to  zero.  This  assumption  makes 
the  calculation  simpler,  although  it  would  make  no 
difference  in  the  case  of  the  Bernoull i-Euler  beam  where  the 
bending  slope  and  the  total  slope  are  the  same. 

The  results  of  the  harmonic  response  problems  are 
presented  in  the  form  of  resonance  curves.  Infinite 
boundary  conditions  are  also  considered  together  with  the 
ordinary  boundary  conditions  such  as  hinged  or  clamped. 

In  order  to  show  the  effect  of  rotary  inertia  and  shear 
deformation,  the  results  are  obtained  for  the  three 
slenderness  ratios;  r=0,  0.04  and  0.08.  When  this 
slenderness  ratio  is  set  to  zero  the  cross-sectional 
dimension  is  in  fact  neglected  and  hence  the  results  become 
those  for  the  Bernoul 1 i -Euler  beam. 

Little  work  has  been  done  for  the  Timoshenko  beam  on 
multiple  supports,  mainly  due  to  the  complexity  of  the 
problem.  The  analytical  solution  by  the  present  method 
therefore  can  not  be  confirmed  by  other  published  results. 
However,  the  results  for  some  extreme  cases  of  zero 
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slenderness  ratio  or  the  single  span  beam  may  be  compared 
with  those  in  [2],  [5],  [9]  and  [10].  For  this  purpose  the 
spring  constants  and  other  parameters  are  chosen  as  those 
found  in  these  references. 


B.  Analytical  Formulation  of  Transfer  Matrix 


Timoshenko  beam  differential  equation 

Consider  a  continuous  beam  shown  in  Figure  4.1(b).  The 
elastic  constraints  provided  at  each  support  are 
characterized  by  a  translational  spring,  kt  and  a  rotational 
spring,  kr  .  All  the  interior  supports  are  assumed  to  be 
identical  maintaining  the  periodicity  of  the  structure. 


(a)  Knife-edge  supports. 
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Figure  4.1.  Uniform  beams  resting  on  equi-spaced 
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The  coupled  equations  for  the  total  deflection,  y  and 
\p ,  the  slope  due  to  bending  within  a  span  are  due  to 
Timoshenko[ 2  1  ] 


EIP+k'GA(  =  0 


_k » GA (  _  M  )_ 

pA at2  k  ax2  ax 


o 


(4.1) 

(4.2) 


in  which 

A  =  cross-sectional  area 
E  =  Young’s  modulus  of  elasticity 
G  =  shear  modulus 

I  =  area  moment  of  inertia  of  cross-section 
k' =  shape  factor  for  cross-section 
p  =  mass  density  of  material. 

Setting 


y  =  7Yexp(icjt) 

(4.3) 

=  4,exp(  icjt. ) 

(4.4) 

{  =  x/7 

(4.5) 

where 

Y  =  normal  function  of  y 
$  =  normal  function  of  ^ 

£  =  nondimensional  local  coordinate 
cj  =  angular  frequency 

=  length  of  a  beam  segment, 
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yields 


s2«P"-(  1-b2r  2s2  )«P+Y'  =  0  (4.6) 

Y"+b2s2Y-¥'=  0.  (4.7) 

Eliminating  $  or  Y  from  equations  (4.6)  and  (4.7),  the 
following  differential  equations  are  obtained 

Y>v +b2 (r2+s2 )Y"-b2 ( 1-b2r2s2 )Y=  0  (4.8) 

4,IV  +  b 2  (  r  2  +  s  2  )4J,'-b2  (  1-b2r2s2  )  4*=  0.  (4.9) 


In  equations  (4.6)-(4.9)  the  nondimensional  parameters  b,  S, 
r  are  defined  as 


b2  = 


(4.10) 


r 2  = 


I 

A/2 


(4.11) 


s  2 


E  I 
k  ’  G  A  7 2 


(4.12) 


and  the  prime  indicates  the  differentiation  with  respect  to 
the  nondimensional  local  coordinate  £. 

Following  Huang's  procedure  further,  the  normal  mode 
solutions  to  equations  (4.8)— (4.9)  are  found  as 


Y  =C 1coshp£+C2sinhp£+C3cosq£+CaSinq£ 


(4.13) 
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$  =D ! coshp^  +D2  s i nhp£+D3cosq£+D4  s inq£ 


(4.14) 


where 


P 

Q 


(4.15) 


The  argument  of  the  hyperbolic  normal  functions,  p  in 
equation  (4.15)  can  become  imaginary  number  when  the  term  in 
the  square  root  bracket  becomes  negative.  This  extension  is 
applied  to  all  the  consequent  equations  involving  with  p. 

The  coefficients  Clf  C2,..  and  Dw  D2  ,  .  .  are  related  by 
equations  (4.6)  or  (4.7)  as 


D  i  —  Q.  C  2 


D 2  —  ttC  i 


(4.16) 


D3  =  j8C* 


Du  =  3 


where  a  and  are  the  nondimensional  constants  defined  by 
a  =  p+b 2  s  2 /p ,  £=q-b2s2/q. 

From  the  beam  flexure  theory  the  nondimensional  shear 
force  V  and  the  bending  moment  M  can  be  expressed  as  [21] 


V  =  -  (k,GA(9y/9x-^) } / 2 /El 


(4.17) 


M 


/  d\p/dx . 


(4.18) 


Substituting  equations  (4. 1 )— (4.5)  into  equations 
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(  4  .  1  7  )  -  (  4  .  1  8  )  the  generalized  forces  are  obtained  as 


V  =  ( b 2  r  2 )  exp  (  icjt ) 
M  =  exp(  io)t )  . 


(4.19) 


Having  the  expressions  for  the  generalized 
displacements,  equations  ( 4 . 1 3 ) - ( 4 . 1 4 )  ,  and  forces, 
equations  ( 4 . 1 9 ) - ( 4 . 20 ) ,  it  is  possible  to  define  a  state 
vector  as 


\  Y  ' 

Y 

=  ( 

* 

* 

V 

b2  r  24»+*i'" 

.  M. 

(4.20) 


Here  the  common  factor  exp(io)t)  was  omitted,  and  all  the 
elements  of  the  vector  {v}  are  nondimensional  quantities. 
Using  the  normal  mode  solution,  equation  (4.13)  and  (4.14), 
the  generalized  displacements  and  forces  at  the  left  and 
right  ends  of  the  beam  segment  are  expressed  in  terms  of 
unknown  constants  {C}. 


■ Y(i-O) ' 

-  1  0  1  0  • 

Cl 

+(£=0) 

.  _ 

’ 

0  a  0  ft 

C  2 

Y(|  =  1  ) 

Ch  Sh  C  S 

1 

C  3 

[  +( {= 1 ) . 

.  aSh  aCh  -ySS  ySC  - 

1  C,  J 

(4.21  ) 


' V(£=0) ’ 

-  0 

qa^ 

0 

-payS 

M(£=0) 

►  ” 

-pa 

0 

0 

V(£  =  1  ) 

qaySSh 

qaySCh 

pa/JS 

-pa/SC 

1  m(£=i  ) 

_  paCh 

paSh 

-qflC 

-qySS  . 

(4.22) 


where  Ch,  Sh,  C  and  S  denote  the  four  normal  functions, 
coshp,  sinhp,  cosq  and  sinq,  respectively,  and  the 
relationship 
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p2+b2r2  =  qp 
q2-b2r2  =  pa 


(4.23) 


was  used  to  simplify  the  expression. 

Impedance  matrix 

The  impedance  matrix  which  relates  the  generalized 
displacements  and  forces  is  then  obtained  by  substituting 
for  the  constants  { C }  in  equation  (4.22)  with  the 
displacements  obtained  from  equation  (4.21)  as 


[Z]  = 


.g-.a.f.g 


(a2  -  f$2  )ShS  +  2a£(  1-ChC) 


"z  ,  , 

Z  1  2 

Z  1  3 

Z  i  «  " 

(4. 

z  2  1 

Z  2  2 

Z  2  3 

Z  2  « 

Z  3  1 

Z  3  2 

Z  3  3 

z  3  « 

-  Z  4  i 

Z  2 

Z  u  3 

Z  «  ti  . 

where  z  t  i  denote  the  components  of  [Z] 


Zyy  =  afi ( aShC+fiChS ) 

z  y  2  =  a^S{ShS+(q/S-pa)  (  1-ChC)  }/(p2+q2  ) 
z  y 3  =  -afi ( aSh+fiS ) 
z  y  it  =  aP  ( Ch~C  ) 

z  2  2  =  aChS-pShC  (4.25) 

z  2  3  =  ( Ch“C ) 

z  2  tt  =  jSSh-aS 


z  3  3 


z  1  1 
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Z3  «  ~  zl  2 
Z  «  4  =  Z  2  2  • 


The  same  impedance  matrix  would  be  obtained  by  the  standard 
finite  element  formulation  procedure,  with  the  interpolation 
functions  being  taken  as  the  four  exponential  functions  or 
two  trigonometric  and  two  hyperbolic  functions.  However, 
the  method  described  in  this  section  is  much  simpler  in  the 
case  of  an  analytical  impedance  matrix. 

The  forces  due  to  springs  (translational  or  rotational) 
are  not  included  in  equation  (4.25),  however  the  effect  of 
springs  at  each  end  of  the  subsystem  can  be  included  easily 
by  adding  to  equation  (4.25)  the  matrix 


[Z]  = 


(4.26) 


where  k2  are  nondimens ional  spring  constants  defined  as 


k,  =  k  t  7 3 /2EI , 
k2  =  kr7/2EI. 


(4.27) 


Transfer  matrix  and  its  inverse 

Once  the  impedance  matrix  is  obtained,  the  expression 
for  the  transfer  matrix  can  be  easily  derived  by  the  simple 
algebraic  procedure  described  in  the  previous  chapter.  The 
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final  expression  of  the  transfer  matrix  for  the  Timoshenko 
beam  element  shown  in  Figure  4.1(b)  is  obtained  as 


[T] 


"tii 

t  1  2 

t  1  3 

t  1  4 

1 

t2, 

t  2  2 

t  2  3 

t  2  4 

P2+q2 

t3 1 

t32 

t  3  3 

t  3  4 

Lt,  i 

t  4  2 

t  4  3 

t  4  4 

where 

tii  =  q)SCh+paC  -  k,(Sh/a-S/£ 
t 1 2  =  pSh+qS  +  k2 (Ch-C) 
t  ,  3  =  (0Sh-aS)/a0 
t , 4  =  Ch-C 

t2  i  =  a£ (qSh-pS )  -  k,  (Ch-C) 
t22  =  paCh+q/?C  +  k2(aSh+0S) 

1 2  3  =  Ch-C 

t 2 u  =  aSh+/iS 

(4.28) 

t3i  =  a0(q2£Sh+p2aS )  -  2k , (q£Ch+paC ) +k ,  2 ( Sh/a-S/0) 
t32  =  b2 (Ch-C )  -  k, (pSh+qS)+k2a0(qSh-pS)-k,k2(Ch-C) 
t33  =  q£Ch+paC  -  k ,  ( Sh/a-S/fi) 

1 3  4  =  a/?  (qSh-pS  )  +  k,(Ch-C) 

1 4  i  =  b2  (Ch-C  )  +  k  !  ( pSh+qS  )-k2  a/?  (qSh-pS)+k  ,  k2  (Ch-C) 
ta2  =  p2  aSh-q 2  /?S  -  2k  2  ( paCh+q0C ) -k  2  2  ( aSh+0S  ) 

1 4  3  =  pSh+qS  -  k2(Ch-C) 

1 4  4  =  paCh+q/?C  +  k  2  ( aSh+/JS  )  . 


This  transfer  matrix  can  also  be  obtained  directly  from  the 
normal  mode  functions  and  beam  flexure  theory  without 
formulating  the  impedance  matrix  as  shown  in  Appendix  B. 
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C.  Special  Cases 

As  mentioned  previously,  the  general  transfer  matrix 
for  the  Timoshenko  beam,  equation  (4.28),  can  be  reduced  to 
that  for  Bernoull i-Euler  beam  by  simply  putting  the 
slenderness  ratio  to  zero,  i.e.,  r2=s2=0,  so  that 


p=q=a=/J=/b . 


(4.29) 


The  transfer  matrix  for  the  Bernoull i-Euler  beam  on 
spring  supports  then  becomes  a  4x4  matrix: 


[T] 


’tii 

1 1 

2 

t  1 

3 

t, 

4 

1 

t  2  1 

t  2 

2 

t  2 

3 

t  2 

4 

2 

t  3  1 

t  3 

2 

t  3 

3 

t  3 

4 

Ltd  i 

t  4 

2 

t  4 

3 

t  4 

4  . 

where 

1 1  i 
1 1  2 
t  1  3 

t  1  4 

t  2  1 
t  2  2 

t  2  3 
t  2  4 
t3i 

t  3  2 
t  3  3 
t  3  4 
t  4  1 


p3 ( Ch+C )  ~ 
p2 ( Sh+S  )  + 
Sh-S 
p(Ch-C) 
p2 ( Sh~S  )  - 
p3 (Ch+C) 
p(Ch-C) 
p2 (Sh+S) 
p‘ ( Sh+S )  - 
p5 ( Ch-C )  - 
p3 (Ch+C )  - 
p4 ( Sh-S )  - 
p5 (Ch-C )  - 


k ! (Sh-S) 
k2p(Ch-C) 


k,p(Ch-C) 


(4.30) 

2k1p3(Ch+C)+k12(Sh-S) 

k,p2 ( Sh+S )+k2p4( Sh-S )-k1k2p(Ch-C) 

k , (Sh-S) 

k1p(Ch-C) 

k 1 p2 ( Sh+S ) +k 2p4 (Sh-S) -k , k2p(Ch-C) 
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t  u  2  =  p 4 ( Sh-S )  +  2k 2p3 (Ch+C ) +k 2  2p2 ( Sh+S ) 
ta3  =  P2 (Sh+S)  +  k2p(Ch-C) 

1 4  4  =  p3 (Ch+C)  +  k  2pz ( Sh+S ) . 


A  special  case  of  this  development  is  the  beam  on 
knife-edge  supports,  which  was  considered  by  Miles[16]  and 
Ayre  and  Jacobsen[1].  The  transfer  matrix  in  this  case  can 
be  obtained  by  setting  the  total  deflection  y  to  zero  on  the 
supports.  The  transfer  matrix  obtained  from  equations 
(4.28)  and  (4.30)  by  this  operation  as  shown  in  Appendix  C 
is  for  the  Timoshenko  beam 


jSSh-aS  * 

‘^ShC-aChS-kzA  A 

_(p2+q2 )ShS-2k2 ( £ShC-aChS ) +k 2  2 A  0ShC-aChS- k 2 A  . 


(4.31  ) 


where  A= { ( a  2 -0 2 ) ShS+2a£ ( 1 -ChC ) } / ( p2 +q 2 ) . 
For  the  Bernoul 1 i -Euler  beam 


(4.32) 

'ShC-ChS-k2 ( 1-ChC)/p  .  (1-ChC)/p 

-2pShS-2k 2 ( ShC-ChS )  +  k 2  2 ( 1-ChC)/p  ShC-ChS-k2 ( 1-ChC)/p  J  . 


Using  the  equations  ( 2 . 25  ) - ( 2 . 26 )  the  total  transfer 
matrices  for  the  N~ span  beams  can  be  obtained  in  terms  of 
the  propagation  constant  and  the  normal  mode  functions 
either  for  the  Timoshenko  or  the  Bernoull i-Euler  beams. 
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D.  Natural  Frequencies 

Timoshenko  beams  on  knife-edge  supports 

Figure  4.2  shows  a  uniform  beam  resting  on  the 
knife-edge  supports  with  a  rotational  spring  provided  at 
each  support.  The  propagation  constants  are  obtained  for  a 
beam  element  with  three  different  k2  values;  A2=0,  1  and 
100,  and  shown  in  Figures  4.3-5.  From  these  propagation 
constant  curves  and  by  the  graphical  techniques  as  explained 
in  the  section  III.D,  the  natural  frequencies  are  obtained 
for  two  different  boundary  conditions;  hinged  and  clamped. 
The  first  four  natural  frequencies  for  the  four-span  beams 
are  listed  in  Table  4.1.  The  corresponding  mode  shapes  are 
shown  for  the  hinged  beam  case  in  Figure  4.6.  Shown  is  the 
total  deflection  as  well  as  that  due  to  bending  alone. 


(a)  Hinged-hinged  boundary. 


2k: 


2k&  2kg  k<g 


"2T 


(b) 


Clamped-clamped  boundary. 


Figure  4.2.  Free  vibrations  of  the  uniform  beams  on 


knife-edge  supports. 
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FREQUENCY  i/b 


Figure  4.3.  Propagation  constant  with  k2=0 
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FREQUENCY  /b 


Figure  4.4.  Propagation  constant  with  k2=1 
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FREQUENCY  /b 


Figure  4.5.  Propagation  constant  with  k 2  =  1 0 0 
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Table  4.1.  Natural  frequencies  of  a  four  span  beam  on 

knife-edge  supports,  /b=  { pA 7  4cj2 /El  }V 


(a)  k2  =  0. 


r  =  0 

0.04 

0.08 

MODE 

HINGE 

CLAMP 

HINGE 

CLAMP 

HINGE 

CLAMP 

1 

3.142 

3.393 

3.083 

3.301 

2.940 

3.090 

2 

3.393 

3.926 

3.301 

3.759 

3.090 

3.412 

3 

3.926 

4.464 

3.759 

4.199 

3.412 

3.709 

4 

4.464 

4.730 

4.199 

4.402 

3.709 

3.833 

(b)  k  2  =  1 . 


II 

O 

0.04 

0.08 

MODE 

HINGE 

CLAMP 

HINGE 

CLAMP 

HINGE 

CLAMP 

1 

3.397 

3.596 

3.316 

3.486 

3.126 

3.241 

2 

3.596 

4.042 

3.486 

3.860 

3.241 

3.492 

3 

4.042 

4.505 

3.860 

4.233 

3.492 

3.733 

4 

4.505 

4.730 

4.233 

4.402 

3.733 

3.833 

( c )  k  2  =  100. 


o 

ii 

0.04 

0.08 

MODE 

HINGE 

CLAMP 

HINGE 

CLAMP 

HINGE 

CLAMP 

1 

4.642 

4.654 

4.337 

4.347 

3.799 

3.804 

2 

4.654 

4.685 

4.347 

4.370 

3.804 

3.816 

3 

4.685 

4.717 

4.370 

4.392 

3.816 

3.828 

4 

4.717 

4.730 

4.392 

4.402 

3.828 

3.833 
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MODE 


2. 


3. 


4. 


TOTAL  DEFLECTION 
BENDING  DEFLECTION 


Figure  4.6.  Normal  modes  of  the  four  span  Timoshenko  beam 


on  knife-edge  supports. 
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Timoshenko  beams  on  elastic  springs 

Figure  4.7  shows  the  beam  supported  on  elastic 
translational  and  rotational  springs.  Natural  frequencies 
are  found  for  three  different  ky  values;  k^=0,  0.5  and  50, 
and  for  two  different  boundary  conditions,  hinged  and 
clamped.  The  rotational  spring,  k2  is  assumed  zero  for  all 
the  results  presented,  although  it  can  be  easily  included. 
Again  the  first  four  natural  frequencies  given  in  Table  4.2 
are  for  the  four-span  beams.  In  Figure  4.8  the 
corresponding  mode  shapes  are  shown  for  the  hinged  beam 
case.  Shown  is  the  total  deflection  as  well  as  the  portion 
of  the  deflection  due  to  bending  alone. 


(a)  Hinged-hinged  boundary. 


(b)  Clamped- clamped  boundary. 


Figure  4.7. 

supports . 


Free  vibration  of  the  beams  on  spring 


. 
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Table  4.2.  Natural  frequencies  of  a  four  span  beam  on 
spring  supports,  /b=  { pA 7  *a>2/EI  }1/4. 


(a )  k  ,  =  0  . 


II 

O 

0.04 

0.08 

MODE 

HINGE 

CLAMP 

HINGE 

CLAMP 

HINGE 

CLAMP 

1 

0.785 

1.183 

0.784 

1.176 

0.782 

1.159 

2 

1.571 

1.963 

1  .563 

1.940 

1  .542 

1  .879 

3 

2.356 

2.749 

2.331 

2.694 

2.264 

2.558 

4 

3.  142 

3.534 

3.085 

3.446 

2.94  1 

3.185 

(b)  k  , 

=  0.5. 

r  =  0 

0.04 

0.08 

MODE 

HINGE 

CLAMP 

HINGE 

CLAMP 

HINGE 

CLAMP 

1 

1 .084 

1.311 

1  .083 

1.307 

1.081 

1 .294 

2 

1.631 

1 .996 

1.623 

1.974 

1.605 

1.917 

3 

2.373 

2.766 

2.349 

2.711 

2.284 

2.577 

4 

3.141 

3.518 

3.082 

3.438 

2.941 

3.185 

(c )  k  ! 

=  50. 

r  =  0 

0.04 

0.08 

MODE 

HINGE 

CLAMP 

HINGE 

CLAMP 

HINGE 

CLAMP 

1 

3.028 

3.045 

2.985 

3.005 

2.872 

2.898 

2 

3.042 

3.065 

3.003 

3.019 

2.897 

2.898 

3 

3.075 

3.608 

3.025 

3.507 

2.898 

3.288 

4 

3.142 

3.689 

3.083 

3.562 

2.940 

3.314 
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MODE 


Figure  4.8.  Normal  modes  of  the  four  span  Timoshenko  beam 


on  spring  supports. 
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E.  Steady-state  Harmonic  Responses 

Timoshenko  beams  on  knife-edge  supports 

A  four-span  beam  with  both  ends  hinged  and  the  beam 
with  infinitely  many  spans  are  shown  in  Figure  4.9.  The 
bending  moment  is  introduced  at  the  center  node  and  the 
bending  slope  at  the  same  point  is  found.  The  results  are 
then  shown  by  the  resonance  curves  -  Figures  4.10  and  4.11. 


(a)  Four-span  hinged-hinged  beam. 


(b)  Infinite-infinite  beam. 


Figure  4.9.  Forced  vibration  of  the  uniform  beams  on 


knife-edge  supports. 
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Figure  4.10.  Resonance  curves  of  four-span  hinged  beams, 
(a)  k 2  = 


0 
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Figure  4.10.  continued... 
(b)  k 2  = 


1 


, 
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Figure  4.10.  continued... 


(c )  k2 


100 
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Figure  4.11.  Resonance  curves  of  infinite  beams  on 
knife-edge  supports. 

(a)  k 2  = 


0 


' 


81 


^1 X 


Figure  4.11.  continued... 


(b)  k 2 


1 
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Figure  4.11.  continued... 


(c )  k2 


100 


■ 
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Timoshenko  beams  on  elastic  springs 

Harmonic  responses  are  obtained  for  the  case  of 
vertical  force  input  and  deflection  output.  The  results  are 
then  shown  by  the  resonance  curves  -  Figures  4.13  and  4.14. 

A  four-span  beam  with  both  ends  hinged  and  an  infinite  beam 
are  taken  assuming  that  the  excitation  is  at  the  center  node 
of  the  beams  as  shown  in  Figure  4.12. 


(a)  Four-span  hinged-hinged  beam. 


|e!wt 


(b)  I nf in i  te- inf  ini te  beam. 


Figure  4.12.  Forced  vibration  of  the  uniform  beams  on 
spring  supports. 
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FREQUENCY  £ 


Figure  4.13.  Resonance  curves  for  a  four  span 
hinged-hinged  beams  on  spring  supports. 

( a )  k  i  =  0  • 
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Figure  4.13.  continued... 


(b)  k ! 


0.5 


RECEPTANCE 
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Figure  4.14.  Resonance  curves  for  an  infinite  beam  on 
spring  supports. 

(a )  k  !  = 


0 
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Figure  4.14.  continued... 


(b)  k  , 


0.5 
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F.  Numerical  Results  by  Finite  Element  Method 

In  order  to  compare  the  results  analytically  obtained 
with  those  from  the  Finite  Element  formulation,  a  span  of 
the  Timoshenko  beam  on  intermediate  supports  is  subdivided 
into  a  number  of  elements.  The  beam  may  thus  be  considered 
an  example  of  the  compound  periodic  structure  -  see  Figure 
4.15. 


(a)  Element  division. 


(b)  Craggs'  element  showing  the  nodal  displacements. 


Figure  4.15.  Timoshenko  beam  element  for  the  Finite  Element 
Method . 


Craggs’  Timoshenko  beam  element  [4]  uses  four  nodal 
degrees  of  freedom,  namely,  the  total  deflection,  Y,  bending 
slope,  V  and  their  first  derivatives  at  each  end  of  the 
element.  In  between  the  nodes  the  total  deflection  and  the 
bending  slope  are  interpolated  by  a  cubic  function  from  the 
nodal  values.  Through  the  general  Finite  Element  procedure 
the  impedance  matrix  is  then  obtained  for  the  beam  element 
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in  the  nondimensional  form  as 


[Z]  =  [K]  -  b 

2  [M] 

(4 

where  the  stiffness 

matrix,  [K]= 

1 

420s2 

'504  42  210 

42  -504 

42 

210 

-42 

56  -42 

0  -42 

-14 

42 

-7 

504s  2  + 1  56 

42s  2  +22  -210 

42 

-504s  2  +  54 

42s  2  -  1 3 

56s  2  +4  -42 

7 

-42s  2  +  1  3 

-  1  4s  2 -3 

504 

-42 

-210 

42 

56 

-42 

0 

Symmetric . 

504s  2  +  1 56 

-42s  2 -22 

56s  2  +4 


and  the  mass  matrix, 


[M]  = 


420 


156 


22 

0 

0 

54 

-13 

0 

0 

4 

0 

0 

13 

-3 

0 

0 

1  56r  2 

22r  2 

0 

0 

54  r  2 

- 1  3r 

4r  2 

0 

0 

1  3r  2 

-  3  r 

156 

-22 

0 

0 

4 

0 

0 

Symmetric 


1  56r  2  -22r  2 
4r  2 


33) 


Here  the  nondimensional  parameters 


are  as  defined  in 
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equations  (  4  .  1  0  )  -  ( 4  .  1  2  )  . 

Unlike  the  analytically  derived  impedance  matrix  in 
equation  (4.24),  this  expression  for  the  impedance  matrix  of 
size  8x8  can  not  be  used  to  represent  the  Bernoul 1 i -Euler 
beam  element  by  setting  r=s=0.  Impedance  matrix  for  the 
Bernoull i-Euler  beam  is  of  size  4x4  and  may  be  found  in  most 
texts  of  Finite  Element  Method. 

At  a  given  frequency,  b,  this  impedance  matrix  can  be 
evaluated  numerically  and  the  transfer  matrix  can  be  formed 
using  equation  (2.8).  This  element  transfer  matrix  is  then 
multiplied  to  form  the  transfer  matrix  relating  the  state 
vectors  at  the  supports.  If  the  beam  is  supported  on  the 
knife-edges  as  in  the  case  of  Figure  4.15,  the  transfer 
matrix  for  a  span  can  be  reduced  in  size  by  the  procedure 
explained  in  the  Appendix  C.  From  the  reduced  transfer 
matrix  thus  obtained  the  eigenvalues  are  calculated  and  the 
logarithmns  of  these  eigenvalues  yield  the  propagation 
constants.  In  Figure  4.16  the  propagation  constant  for  the 
flexural  motion  is  plotted  against  frequency  for  the  two 
slenderness  ratios,  r=0.04  and  0.08.  For  each  slenderness 
ratio  a  span  of  the  beam  is  divided  into  one,  two  and  three 
elements.  From  these  propagation  constant  curves  the 
natural  frequencies  are  obtained  by  the  graphical  method  as 
explained  in  Section  III.D.  Table  4.3  gives  the  results  for 
the  hinged-hinged  beams  in  terms  of  relative  errors  compared 
with  the  exact  solutions. 


PROPAGATION  CONSTANT 
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Figure  4.16.  Propagation  constant  by  the  Finite  Element 
Method.  (a)  r=0.04. 
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Figure  4.16.  continued... 


(b)  r=0 . 08 . 
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Table  4.3.  Errors  of  Finite  Element  results  for  the 

four  span  Timoshenko  beam  on  knife-edge  supports. 


(a )  r  = 

• 

o 

• 

o 

(  %  ) 

NUMBER  OF 

ELEMENTS 

PER  SPAN 

GROUP 

MODE 

/b 

1 

2 

3 

1 

3.083 

4.7 

o 

• 

o 

o 

• 

o 

1 

2 

3.301 

8.7 

0.3 

0.  1 

3 

3.759 

25.3 

0.7 

0.3 

4 

4.199 

64.9 

1  .5 

0.7 

5 

5.881 

— 

1  .5 

0.  1 

2 

6 

6.046 

- 

2.6 

ID 

• 

O 

7 

6.378 

- 

5.2 

1  .2 

8 

6.689 

- 

9.4 

2.  1 

(b)  r  = 

o 

• 

o 

00 

• 

(  %  ) 

NUMBER  OF 

ELEMENTS 

PER  SPAN 

GROUP 

MODE 

/b 

1 

2 

3 

1 

2.940 

1.5 

o 

• 

o 

o 

• 

o 

1 

2 

3.090 

4.0 

0.5 

0.3 

3 

3.412 

13.6 

1  .6 

1  .0 

4 

3.709 

38.5 

2.9 

1  .7 

5 

5.192 

- 

0.6 

0.1 

2 

6 

5.260 

- 

1.9 

0.7 

> 

7 

5.400 

- 

4.9 

o 

• 

eg 

8 

5.520 

- 

8.6 

3.5 
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G.  Discussion 

It  can  be  easily  shown  that  the  same  impedance  matrix 
as  in  equation  (4.7)  can  be  obtained  by  the  classical  Finite 
Element  Method  using  the  four  normal  modes  as  the 
interpolation  functions,  instead  of  the  polynomials. 

However,  this  procedure  requires  more  calculation  than  that 
is  needed  by  the  unknown  coefficient  method  employed  in  this 
chapter.  If  one  requires  an  approximate  solution,  one  can 
use  the  polynomials  (most  commonly  cubics)  to  obtain  the 
impedance  matrix  in  the  familiar  form  as  in  equation  (2.2). 
The  precise  method  of  formulating  the  impedance  matrix  is 
not  important  in  the  present  study  because  the  present  study 
starts  from  the  impedance  matrix.  One  should  be  cautious  of 
the  sign  convention  used  in  the  present  study  because  it  is 
consistent  with  that  of  Finite  Element  Method  not  that 
commonly  used  in  many  texts  -  see  equation  (4.18)  for 
example . 

L i n [  10]  used  considerable  effort  to  find  the  expression 
for  [Tn]  for  the  Bernoulli-Euler  beam  on  spring  supports. 
While  this  approach  can  also  be  used  for  the  Timoshenko  beam 
case  the  Cay ley-Hami 1 ton  theorem  and  the  similarity 
transform  methods,  shown  are  more  efficient.  Besides,  the 
present  method  is  aimed  at  the  generally  shaped  periodic 
structures  not  just  uniform  beams. 

As  mentioned  previously,  the  four  components  of  the 
transfer  matrix  for  the  mono-coupled  systems  are  in  fact  the 
frequency  equations  for  a  single  span  system  with  different 
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boundary  conditions.  Therefore  three  frequency  equations 
can  be  obtained  for  the  Timoshenko  beams  with  three 
different  boundary  conditions  from  the  transfer  matrix  in 
equation  (4.31);  tii  (or  t22)  for  clamped-hinged ,  t12  for 
c lamped-c lamped  and  t21  for  hinged-hinged  case.  These 
frequency  equations  are  the  same  as  those  found  in  [7] 
except  in  the  case  of  hinged-hinged  boundary,  where  the  two 
equations  yield  the  same  set  of  solutions  in  the  low 
frequency  range.  However,  in  the  high  frequency  range  where 
the  argument  of  the  hyperbolic  normal  functions,  p  becomes 
imaginary  number  so  that  these  functions  in  fact  become 
trigonometric,  Huang’s  equation  fails  to  give  all  the 
frequency  roots. 

All  the  results  are  obtained  for  the  three  slenderness 
ratios;  r=0,  0.04  and  0.08.  The  nondimensional  parameter,  s 
is  set  to  twice  these  values  assuming  that  the  shape  factor, 
k'=2/3  (or  0.65)  and  E/G=8/3  (or  2.6),  which  was  used  in  the 
references . 

Figures  4. 3-4. 5  show  the  propagation  constants  for  the 
Timoshenko  beam  on  the  knife-edge  supports  and  with  the 
rotational  springs,  k2=0,  1  and  100.  With  r=0  and  k2=0  the 
propagation  constant  curve  becomes  identical  to  Figure  3.5 
and  also  that  in  [20]  for  the  Bernoull i-Euler  beam.  The 
propagation  constant  has  a  repeating  regular  pattern  of 
interval  2tt,  with  the  propagation  zones  lying  in  the 
frequency  ranges  m7r</b<  (m+0 . 5  )  tt  ,  approximately.  Therefore, 
once  the  natural  frequencies  of  the  first  group  are 
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obtained,  the  higher  group  natural  frequencies  are  simply 
obtained  by  adding  multiples  of  7 r  to  these  frequencies. 

This  regular  pattern  becomes  distorted  as  the  slenderness 
ratio  increases  by  shifting  the  curve  towards  the  lower 
frequency  side,  that  is,  the  natural  frequency  becomes  lower 
with  increasing  slenderness  ratio.  Above  a  certain 
frequency  the  alternating  pattern  of  attenuation  and 
propagation  zones  completely  breaks  and  there  can  be  an 
additional  natural  frequency  in  a  group. 

The  increase  in  the  rotational  spring  constant,  k2 
makes  the  propagation  zone  narrower  by  moving  the  lower 
bounding  frequency  towards  the  upper  bounding  frequency 
giving  higher  natural  frequencies.  Since  the  upper  bounding 
frequencies  are  the  natural  frequencies  of  the  single  span 
clamped-clamped  beam,  these  are  independent  of  the 
rotational  spring,  k2.  The  lower  bounding  frequencies,  on 
the  other  hand,  which  are  the  natural  frequencies  of  the 
single  span  hinge-hinged  beam  for  which  the  ends  can  rotate, 
increase  for  the  higher  k2  .  As  k2  increases  further  these 
natural  frequencies  approach  the  upper  bounding  frequencies 
of  the  propagation  zones  which  are  the  natural  frequencies 
of  the  single  span  beam  with  both  ends  clamped.  The 
propagation  of  the  flexural  waves  is  then  possible  only  at 
these  frequencies  and  every  span  of  the  beam  becomes 
virtually  isolated. 

In  the  case  of  beams  on  spring  supports,  little  is  to 
be  gained  by  considering  the  propagation  constant  and  the 
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natural  frequencies  are  found  by  iteratively  solving  the 
frequency  equations.  The  results  are  shown  in  Table  4.2. 
With  both  k]  and  k2  being  zero  the  beam  becomes  a  single 
span  beam  with  no  intermediate  supports,  which  is  identical 
to  the  beam  considered  in  [7],  The  propagation  constants  in 
this  case  are  M  =  ±p,  ±qi  and  the  natural  frequencies  are  jr/N 
(j=1,2,...)  for  the  Bernoull i-Euler  beam  with  both  ends 
hinged;  and  approximately  (j  +  0.5)7r/N  (j  =  1,2,...)  for  the 
clamped-clamped  case.  The  reason  that  the  number  of  spans, 

N  appear  in  these  frequencies  is  that  the  frequency 
parameter  b  is  based  on  the  beam  element  length. 

By  increasing  the  translational  spring  constant,  ki, 
the  natural  frequencies  increase  approaching  the  natural 
frequencies  of  the  beams  on  knife-edge  supports.  The  effect 
of  the  rotational  spring  would  be  higher  natural 
frequencies,  however,  it  is  not  considered  in  this  case. 

The  total  deflections  shown  in  the  mode  shapes  in 
Figure  4.6  and  4.8  are  identical  to  those  for  the 
Bernoull i-Euler  beams  except  that  the  total  deflection  of 
the  Timoshenko  beam  is  composed  of  deflections  due  to 
bending  and  shear.  When  the  mode  shapes  are  simple 
sinusoidal  curves  as  the  first  mode  in  Figure  4.6  and  all 
the  modes  in  Figure  4.8,  there  exists  a  constant  which 
governs  the  proportion  of  the  deflections  due  to  bending  and 
shear  throughout  the  beam.  The  proportion  of  the  shear 
deformation  increases  as  the  frequency  increases  and  at  the 
sufficiently  high  frequency  normal  modes  of  pure  shear  can 
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appear.  The  translational  or  the  rotational  springs  at  the 
supports  do  not  affect  the  mode  shapes. 

The  effect  of  the  rotary  inertia  and  shear  deformation 
on  the  natural  frequencies  is  shown  in  the  form  of  a 
correction  factor  given  in  Figure  4.17.  The  corrections  for 
the  single  span  hinged-hinged  beam  are  similar  to  those  for 
clamped-free  beam  as  shown  in  [7],  In  the  case  of  the  four 
span  beam,  the  corrections  for  the  first  four  modes  fall  in 
between  those  of  the  first  and  second  mode  of  the  single 
span  beam. 

There  are  two  resonances  and  two  zero  responses  in  each 
frequency  band  in  Figures  4.10  and  4.11.  A  four-span  beam 
has  four  natural  frequencies,  of  which  two  correspond  to 
symmetric  modes,  the  other  two  to  anti-symmetric  modes. 

Since  the  external  moment  is  applied  at  the  center  node,  it 
only  excites  the  an t i -symmet r ic  modes.  Hence  the  resonance 
occurs  at  these  two  natural  frequencies  corresponding  to 
anti-symmetric  modes  and  the  bending  slope  vanishes  at  the 
center  at  the  natural  frequencies  corresponding  to  symmetric 
modes.  It  is  obvious  then  that  the  number  of  these 
resonances  and  zero  receptances  in  each  group  increases  with 
the  number  of  spans,  making  the  analysis  very  difficult  as  a 
slight  change  in  the  frequency  results  in  the  change  of 
response  from  zero  to  infinite  or  vice  versa.  Therefore  it 
is  more  convenient  to  assume  the  number  of  spans  being 
infinite  in  this  circumstance,  which  not  only  saves  the 
computational  efforts  but  yields  reasonable  results  as  shown 
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(a)  Single  span  beam. 
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Figure  4.17.  Effect  of  rotary  inertia  and  shear 

deformation  on  the  natural  frequencies  of  hinged 
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in  Figure  4.11.  In  the  case  of  infinite  beams  there  is  a 
single  resonance  and  a  zero  receptance  in  each  frequency 
band;  the  resonating  frequency  is  the  natural  frequency  of  a 
single  span  beam  with  both  ends  hinged,  and  zero  receptance 
at  the  natural  frequency  of  a  clamped-clamped  beam  element 
smoothing  the  response  in  between  these  two  bounding 
frequencies.  A  similar  trend  can  be  observed  in  Figures 
4.13  and  4.14,  which  show  the  vertical  deflection  due  to 
vertical  excitation  introduced  at  the  center  nodes  of  a 
four-span  beam  and  an  infinite  beam  on  elastic  springs.  The 
response  does  not  have  to  be  found  only  at  the  point  of 
excitation.  The  state  vectors  at  any  nodes  can  be  obtained 
by  equation  (3.12),  once  the  equation  (3.11)  is  solved, 
since  the  upper  half  of  the  entries  in  the  state  vectors  are 
the  displacement  responses.  The  harmonic  response  problems 
considered  are  for  the  cases  of  force-input 
displacement-output  problems.  However,  since  the  input 
vector  {w}  in  equation  (3.11)  contains  both  the  displacement 
and  force  terms,  the  opposite  case  can  also  be  treated  in 
identical  ways.  In  fact  even  the  mixed  input  problems  can 
be  considered,  however,  it  is  rare  to  find  such  applications 
in  practice. 

It  should  be  noted  that  the  response  at  zero  frequency 
is  in  fact  the  response  due  to  the  static  force.  Hence  the 
present  method  of  analysis  does  not  exclude  the  static 
analysis  of  periodic  structures. 
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Whether  the  impedance  matrix  is  obtained  analytically 
or  by  the  Finite  Element  Method,  the  transfer  matrix 
analysis  is  identically  applied  in  the  present  study.  This 
was  demonstrated  by  using  the  impedance  matrix  in  reference 
[4].  The  accuracy  of  the  Finite  Element  solution  depends 
largely  upon  how  closely  the  interpolation  function  can 
approximate  the  true  shape  of  the  normal  modes.  If  the 
cubic  function  is  used  as  with  the  Craggs'  element,  the 
function  may  not  satisfactorily  approximate  the  deformation 
corresponding  to  the  second  or  higher  group  natural  modes. 
Therefore  only  the  first  group  of  natural  frequencies  are 
obtained  with  a  single  element  per  span.  As  the  number  of 
elements  per  span  increases  the  propagation  constant  curve 
from  the  Finite  Element  more  closely  resembles  that  of  the 
analytically  obtained. 

Generally,  the  accuracy  of  the  Finite  Element  solution 
deteriorates  as  the  frequency  increases.  For  the  multi-span 
beam,  however,  the  accuracy  the  lower  mode  solution  in  the 
next  group  is  better  than  that  of  the  last  mode  in  the 
previous  group.  This  can  be  also  noted  in  the  propagation 
constant  curve.  In  each  propagation  zone  the  discrepancy 
increases  with  the  frequency  but  the  discrepancy  at  the 
lower  bounding  frequency  of  the  next  propagation  zone  is 
smaller  than  that  at  the  upper  bounding  frequency. 

The  same  four  span  Timoshenko  beam  was  also  treated  by 
the  ordinary  Finite  Element  Method.  The  natural  frequencies 
were  obtained  by  finding  the  eigenvalues,  b2  from  the  global 
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impedance  matrix  assembled  from  the  element  impedance 
matrices.  The  solutions  by  this  ordinary  Finite  Element 
Method  are  the  same  as  those  obtained  through  the 
formulation  of  the  transfer  matrix  within  the  round-off 
error.  The  advantage  of  the  present  method  of  analysis  is 
then  the  computational  efficiency.  The  computing  time  (CPU) 
required  for  these  two  methods  is  shown  in  Figure  4.18. 
which  represents  only  the  specific  case  of  the  example  in 
the  present  study.  However,  the  general  trend  can  be 
observed  indicating  that  the  increasing  computational  effort 
is  necessary  to  obtain  better  solutions  by  the  ordinary 
Finite  Element  Method.  On  the  other  hand,  the  computational 
effort  is  virtually  constant  regardless  of  number  of 
elements  per  span  when  the  transfer  matrix  method  is  used. 
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Figure  4.18.  Comparison  of  the  computing  time  used. 


V.  CONCLUDING  REMARKS 


There  have  been  two  different  approaches  used  for  the 
vibration  analysis  of  periodic  structures;  the  matrix  method 
and  the  wave  approach.  The  present  method  is  a  blending  of 
these  two  approaches  putting  more  emphasis  on  the  transfer 
matrix  and  merging  the  properties  of  propagation  constants 
and  the  eigenvalues  of  the  transfer  matrix.  Since  periodic 
structures  can  be  of  complex  shapes  so  that  the  analysis  can 
be  done  by  numerical  means  only,  the  transfer  matrix  in  the 
present  study  is  derived  from  the  impedance  matrix  which 
commonly  formulated  by  the  Finite  Element  Method.  These 
ideas  are  not  really  new  and  had  been  used  by  previous 
authors.  However,  the  present  study  is  somewhat  more 
complete  and  easier  to  implement  since  the  method  suggested 
is  quite  general,  efficient  but  very  simple.  Some  of  the 
advantages  of  the  present  study  which  have  been  discussed  in 
the  previous  chapters  are: 

1.  The  periodic  structure  can  be  of  any  complex  form. 
Many  other  methods  are  only  applicable  to  specific  type 
of  periodic  structure,  but  the  present  study  can  handle 
a  periodic  structure  composed  of  different 
substructures,  or  a  structure  which  is  not  even  periodic 
at  all.  The  periodic  unit  can  be  either  symmetric  about 
its  center  or  unsymmetric. 

2.  The  periodic  structure  may  have  any  type  of  extreme 
boundaries.  Even  the  infinite  structures  can  be  handled 
by  the  present  method,  which  can  not  be  analyzed  by 
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other  matrix  methods.  The  wave  approach  on  the  other 
hand  can  not  satisfactorily  handle  the  ordinary  boundary 
conditions . 

3.  Free  or  forced  vibration  analyses  are  performed  with 
practically  the  same  procedure.  Other  methods  usually 
can  handle  either  free  or  forced  analysis. 

4.  As  many  natural  frequencies  as  needed  or  those  in 
the  specific  frequency  ranges  can  be  obtained  by  the 
present  method.  The  Finite  Element  Method,  on  the  other 
hand,  yields  only  the  same  number  of  natural  frequencies 
as  the  size  of  matrix  starting  always  from  the  lowest 
natural  frequency. 

5.  The  results  can  be  obtained  either  analytically  or 
numerically,  depending  on  what  type  of  impedance  matrix 
used.  Other  studies  are  commonly  directed  to  just  one 
of  these.  The  method  yields  accurate  results  without 
any  complicated  measures  of  avoiding  errors  due  to  the 
digital  computations.  This  is  accomplished  by  reducing 
the  computational  effort  to  a  minimal  amount. 

Any  periodic  structure  can  be  taken  as  the  numerical 
example  for  the  present  study.  However,  it  is  believed  that 
the  results  for  the  Timoshenko  beams  on  multiple  supports 
are  readily  applicable  to  practical  engineering  problems. 
These  results  are  more  useful  for  the  short  and  stubby  beams 
and  for  the  high  frequency  analysis,  where  the 
Bernoulli -Euler  beam  theory  fails  to  give  acceptable 
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In  the  cases  where  the  high  frequency  vibration  is 
concerned  or  the  periodic  structure  contains  many  bays,  the 
damping  effect  may  be  noticeable  and  the  analysis  should 
include  the  damping  terms.  Further  study  may  be  directed  to 
understanding  of  the  wave  phenomena  in  terms  of  the 
propagation  constant  and  the  transfer  matrices.  It  may  be 
also  possible  to  extend  the  application  of  the  present  study 
to  engineering  problems  other  than  the  static  and  dynamic 
analysis  of  the  structures,  such  as  heat  transfer,  fluid 
flow  or  some  other  problems  which  can  be  expressed  in  the 
form  of  difference  equations. 
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VI.  APPENDIX 


A.  Impedance  matrices  for  spring-mass  systems 


1)  System  of  symmetric  elements 

It  can  be  shown  that  the  mass  and  the  stiffness 
matrices  can  be  expressed  for  the  symmetric  element  shown  in 
Figure  3.3(a)  in  the  form 


[M] 


(6.1) 


[K] 


(6.2) 


The  impedance  matrix  is  then  expressed  in  the  nondimens i ona 1 
form  as 


[Z] 


*  1-b2 

-  1 

.  -1 

1-b 

(6.3) 


where  b2=mo>2/k.  Using  this  impedance  matrix  and  the 
procedure  described  in  equation  (2.8)  ,  the  transfer  matrix 
as  in  equation  (3.24)  is  obtained. 

2)  System  of  unsymmetric  elements 

For  the  unsymmetric  element  the  mass  and  the  stiffness 
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matrices  can  be  obtained  as 


[M] 


0  ' 

0  J  , 


(6.4) 


[K] 


(6.5) 


The  impedance  matrix  in  this  case  becomes 


[Z] 


1 -2b2  -1  ‘ 

-1  1  . 


(6.6) 


From  this  impedance  matrix,  the  transfer  matrix  in  equation 
(3.29)  is  obtained. 


1 1 1 


B.  Derivation  of  transfer  matrix  by  direct  method 

The  nondimensional  shear  force  and  the  bending  moment 
shown  in  Figure  4.1  can  be  evaluated  from  equation  (4.19) 
including  the  terms  due  to  springs,  k^  and  k2  as 


V  =  b2  r  24*+4J"  +  k  y  Y 
M  =  -4*'+k24\ 

at  £  =  0  ,  and 


(6.7) 


V  =  b2r24»+^”-k1Y 
M  =  -4s ’  -k 2^ , 


(6.8) 


at  £  =  1  . 

Thus  the  state  vectors  at  £=0  and  £=1  can  be  written  in 
terms  of  the  normal  mode  functions  and  the  coefficients,  {C} 
as 
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The  matrix  in  equation  (6.10)  is  postmultiplied  by  the 
inverse  of  the  matrix  in  equation  (6.9)  to  yield  the 
expression  for  the  transfer  matrix  identical  to  equation 


(4.28)  . 


C.  Derivation  of  the  reduced  transfer  matrix 

The  transfer  matrix  for  the  beam  on  knife-edge  supports 
can  be  obtained  from  that  for  the  beam  on  spring  supports  in 
equation  (4.28)  by  eliminating  the  components  corresponding 
to  the  deflection,  /  and  the  shear  force,  V.  Since  the 
deflection,  Y  is  not  permitted  on  the  knife-edge  supports 


(6.11) 


Y  =  0  =  t12*+t13V+t1ftM, 


where  t jj  denotes  the  /- th  row  and  j-th  column  entry  of  the 
transfer  matrix  in  equation  (4.28)  with  k^=0.  The  third 
column  of  the  transfer  matrix  then  can  be  replaced  by  the 
relationship  obtained  in  equation  (6.11)  leaving  the  reduced 
matrix  of  size  two  as 


[T] 


(6.12) 


The  transfer  matrices  in  equations  (4.31)  and  (4.32)  are 
obtained  by  this  procedure  from  the  transfer  matrices  in 
equations  (4.28)  and  (4.29). 
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